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SUMMARY 


A theoretical and experimental study was made of the buckling 
characteristics of perfect and imperfect circular cylindrical shells subjected 
to dynamic axial loading. 

The experimental tests were performed on a specially designed dynamic 
testing-machine which was capable of producing controlled ramp-type' loads at 
rates' ranging- from the' quasi-static' up to higher than 200,000 pounds /second. 

Ten shell specimens 'were tested including two near-perfect shells , - seven- 
shells' with "axi symmetric sinusoidal imperfections of a variety of ■ amplitudes 
and wave-lengths and one shell with quasi-random axisymmetric imperfections. 

The shells were' produced from a photo-elastic epoxy plastic using a spin-casting 
technique - .' "The 'imperfection profiles were machined into the shell walls using 
a high- precision' hydraulic- tracing apparatus . - For three of 'the shells with 
sinusoidal 'imperfections imperfection profiles were cut • on the inner surface 
alone" while" for "the- 'remaining- four’ shells with sinusoidal imperfections and 
for the" shell 'with quasi-random imperfections , a special manufacturing procedure 
was adopted which produced shells with matching inner and outer profiles, thus 
providing effectively constant thickness walls. 

• Experimental data included dynamic buckling loads (12k data points), 
high speed photographs of the buckling mode shapes and observations of the 
dynamic stability of shells subjected to rapidly applied sub-critical loads. 

A mathematical model was developed to describe the dynamic behaviour 
of perfect and imperfect shells . This model was based on the Donnell-Von K6rm£n 
compatibility and equilibrium equations and had a wall deflection function 
incorporating five separate modes of deflection. Close agreement between 
theory and experiment was found for both dynamic buckling strength and buckling 
mode shapes . 
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NOTATION* 


A,B 


E (w(x,y,t) f (x,y,t) ) 

F g ,F 1’ ’*• F 13 
h 

h 

k 

K,M,N 

k 

cr 


N ,N ,N 
x’ y’ xy 

P 

P 


deflection mode amplitudes 



one-dimensional acoustic wave speed 
(= 6.77 x 10^ inch/ sec for epoxy) 


modulus of elasticity 

equilibrium equation 

AIRY stress function 

coefficients of AIRY stress function 

shell wall thickness 

average shell wall thickness 

constant related to gas flows in loading 
cylinder of testing machine 

axial and circumferential wave-numbers 


12(l-v 2 )' 

L h 2 R 2 J 


ft 


critical axial wave-number 


length of sample employed in split Hopkinson 
bar experiment 

length of shell 

membrane stress resultants 

axial load on shell 

air pressure in loading cylinder of testing 
machine 

initial air pressure in loading cylinder 



2irEh 2 

[3(l-v 2 )] h 
shell radius 
time 


classical buckling load 


u,v,w displacements measured in axial, circumferential 

and radial direction respectively 

V closing velocity of shell end plates 

* Symbols not appearing in this table are defined locally in the text. 


v 




9 * * ' 


radial displacement function amplitude 
coefficients 


W 

o 


W, 


02 


w 


w w 

- 1 ’ ••• —02 


x,y,z 


initial radial imperfection 

imperfection amplitude coefficients 

axial, circumferential and radial coordinates, 
respectively 


Greek Symbols 
a, g 

Y 

6 

x 

A 

e 

X 

T 


y 


V 


p 


a 

x 


I 

X 



radial deflection mode wave-numbers 

non-linear loading rate parameter 

axial compression of shell 

Laplace operator 

strain 

wave-length 

shear stress 

y 

— imperfection amplitude parameter 
h 

Poisson's ratio 
material density 

axial and circumferential stress respectively 
(tensile .= positive) 

axial and circumferential uniform compressive 
stress (compressive = positive) 


Superscripts 

Dot notation (w, etc.) represents differentiation with respect to time. 


Subscripts 

A,B material type symbols (wave transmission) 
cr critical 

cJl classical 

I material interface 

When a subscript follows a comma, this indicates partial differentiation of the 

principle variable with respect to the subscript. 



1 . INTRODUCTION 

■ Thin-walled cylindrical shells are common structural elements . in 
many areas of" engineerings "especially in'the ;aerospace field;' -For the designer 
attempting to build efficient and reliable structures, the stability of these, 
shells under static "and ‘dynamic "loads .is of . extreme importance 

It has been recognized ‘for ■ some time that the buckling loads predicted 
by the "classical ‘small-deflection ‘theory ‘ are ..considerably .higher than - those 
which maybe "attained with ‘all but the most -meticulously prepared experimental 
specimens; - In a ‘study published in ‘1-93^, Donnell-*- ; attempted to resolve , this 
discrepancy "through "the use -of a large-deflection analysis which considered . 
the effects" of ;.initial " imperfections: and plastic material deformations. 

However - ; -this' theory "contained some serious oversimplifications, and the 
results"were "far too unreliable -to be of ..any use to the practicing engineer. 

In 19^1 - Von- Kdrmdn - and - Tsien 2 used . a similar large-deflection theory to show- 
that a stable post-buckled equilibrium state could. exist under axial loads : 
substantially, lower than -the classical buckling load. It was felt that the 
minimum post-buckling load would provide a useful conservative estimate,. of 
the- buckling strength of cylindrical shells and a number of .workers applied 
themselves . to the. task of ref ining the basic Von Kdrmdn-rTsien approach in- 
pursuit of a true lower bound, This search effectively came to, an end in 
1966 when Hoff, Madsen and Mayers® demonstrated .that . for shells with very 
high radius to. thickness ratios, the .minimum post-buckled equilibrium state 
could, exist with essentially zero axial load. 

An alternative approach was 'followed by Koiter, who in, 19*+ 5 
published a general non-linear theory of elastic stability' ^ based on 
consideration ,o-f the effect of initial imperfections on the pre-buckling - 
deformations and, ultimately, on buckling loads. Koiter also prepared a 
special version of the theory 5 which gave a quantitative .relationship .between . 
buckling loads and imperfection geometry for the case of axially loaded, 
cylindrical shells with finite axisymmetr.ic imperfections. This very important - 
work was originally published in Dutch and it was not until 1963 when English 
translations appeared. that the theories became well-known outside of the - 
Netherlands. Since that time, however ,. Koiter ' s work has become widely 
accepted, .both, as a general theory, for the stability of imperfection-sensitive 
structures and -as a useful model in the .study of the buckling of cylindrical 
shells... Several contemporary. workers have generalized and expanded this 
basic approach (References 6 to 10) and experimental evidence has been produced 
which largely confirms the theoretical predications (References 10 to 13). 
Recently a safe design criterion was proposed"* - which took into account the. 
actual imperfections of realistic shell structures. Good reviews of. current 
theory are found, in References 15 , 16 and -IT. 

In addition to the previously ; cited work on static buckling, efforts 
have been made to predict the effect of dynamically applied loads on the 
stability of shells.- Budiansky and Hutchinson^® A9 and Danielson ^ 0 presented 
general theories for .the dynamic buckling of , imperfection-sensitive structures. 
Studies more specifically, related to cylindrical shells have been made by 
Agamirov and Vol'mir^-*-, Coppa and Nash^, Roth and, Klosner^®, Bieniek, Fan 
and Lackman^ , Dym ^ and Dietz^®. These investigations have generally been 
based on some form of the Donnell-Von Karmen large-deflection theory modified 
to include radial inertia. Roth and Klosner . considered a configuration similar 
to that investigated by Budiansky and Hutchinson, that is,, a cylindrical shell 


1 



subjected "to an axial loadrapplied either stepwise or as a pulse of limited ■ 
duration. ' Both "Agamirov and Volhnir aM Bieniek, Fan and Lackiiian studied 
the caae~of ‘ a' dynamically applied external pressure load. Agamirov and Vol'mir 
also considered the problem of a shell loaded axially at a constant rate (ramp 
loading). However, while their paper contains an equation governing the motions 
of the shell wall , no solutions or results were provided; 

- Nash*; '"Dyfti‘ and Dietz have each" investigated the buckling of: shells 

under'a" constant-rate "of ‘axial compression or end shortening. Dynr presented 
a fairly simple - '' analysis involving a single asymmetric deflection mode . - No ' ' ‘ ' 
final ' solution" wa‘S‘ given; rather , the equation of motion was compared to"an‘ " 
earlier 'equation' derived by Hoff "'to describe the 'behaviour 'of an imperfect ' 
column" ‘loaded yin 'a • similar ' fashion. Nash, in the papers co-authored by himself 
and A. P. Coppa, followed a similar approach but based his analysis on a two 
degree ;; of freedom deflection mode of the form: 

W(x,y,t) = A(t)sinoxsin3y + B(t )sin 2 axsin 2 gy . 

The equations of motion of the shell wall derived in this paper were integrated 
using a Runge-Kutta type numerical procedure. By this technique, Coppa was 
able to predict theoretical dynamic buckling loads. Unfortunately, due to 
oversimplification of the mathematical model, especially with respect to the 
assumed deflection modes, the published results are unreliable. A more sophisti- 
cated and generalized form of this theory was given by Dietz who considered 
the case of dynamic buckling of shells with closely spaced eccentric stiffen- 
ing. Here, the constant-thickness shell was considered only as a limiting 
case of the more general configuration and thepublished results were restricted 
to the reinforced cases. A three degree of freedom deflection mode was used 
in this analysis, this being composed of a time dependent term representing 
uniform radial expansion plus two trigonometric terms similar to those employed 
by Coppa. 


Relatively little experimental work has been done on the dynamic 
buckling of cylindrical shells. Agamirov and Vol'mir have mentioned some experi- 
ments which they claim confirmed the basic trends predicted by their theoreti- 
cal work but no specific results or descriptions of these exp er inept s are 
available in the open literature. Coppa^ and Schwieger^l have reported 
experiments in which shell models were buckled under the force of a falling 
weight. In both cases, however, the shells were simply rolled from sheet and 
adhesively joined at the seam, a method which provides very irregular and 
unreliable specimens. Further, in both experiments, the impact loads were 
applied directly to the free ends of the shells. Static theories and tests 
have shown that weak or non-existent end constraints of this nature can 
seriously degrade the load carrying capabilities of the shell, structures. 

In summary, then, although the buckling of circular cylindrical shells 
under static loads is fairly well understood at present, the study of the 
response of these structures to dynamically applied loads is considerably 
less well advanced. Several major deficiencies are apparent. First, the 
theoretical analyses attempted so far have generally made use of rather limited 
imperfection and deflection mode functions and the published results have been 
restricted to a few specific cases so that it is impossible to make general 
conclusions. Further, due to the lack of reliable quantitative experimental 
results, it has not been possible to confirm the theoretical predictions or to 
check the validity of the assumptions impl-ied in the theoretical studies. 
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This report contains the results of an experimental study in which a 
number of carefully manufactured near-perfect and imperfect shell specimens were 
buckled under dynamic axial loads in a specially designed testing machine . It 
also outlines the development of a general mathematical model for dynamic shell 
response and contains a direct comparison between the experimental results and 
the theoretical predictions. Both regular and quasi-random axisymmetric imper- 
fections are considered and experimentally observed buckling modes are compared 
to the optimum deflection functions suggested by the theoretical analysis. The 
mathematical model is based on a modified version of the Donnell-von Karman 
compatibility and equilibrium equations. A set of non-linear equations governing 
the shell wall displacements are developed from these using a Galerkin procedure 
and a solution is found by numerical integration. A five-mode deflection func- 
tion which includes terms representing uniform radial expansion of the shell 
wall, axisymmetric and asymmetric sinusoidal deflections is used. This deflec- 
tion function permits a high degree of flexibility when modelling the behaviour 
of specific shell geometries. 

2. THEORETICAL ANALYSIS 


This section describes the development of a mathematical model for 
prediction of the force-displacement relation and ultimate buckling load of an 
imperfect cylindrical shell subjected to a dynamically applied axial load. 

The development is based on a modified version of the Donnell-von 
Karman compatibility and equilibrium equations. These are essentially a pair 
of fourth- order , non-linear differential equations with time and geometric 
co-ordinates as independent variables. Substitution of a deflection mode 
function and application of a Galerkin procedure transformed these equations 
into a series of coupled, time-dependent differential equations governing the 
radial motions of the shell wall. These equations were then integrated numeri- 
cally to provide a time history of the behaviour of the shell under application 
of a dynamic load. This procedure is similar in concept to several of those 
described in the preceding section. However, the displacement mode function 
applied here was designed to permit a high degree of flexibility in the analysis 
of the effect of specific imperfection geometries. 

2.1 Compatibility and Equilibrium Equations 


The Donnell-von Karman compatibility and equilibrium equations used 
in this analysis were modified to include terms that accounted for the effect of 
radial inertia and the existence of initial imperfections. In this formulation, 
the shell displacements are defined with respect to a Cartesian co-ordinate system 
with its origin on the shel^L.'wall mid-surface, the x and y co-ordinates oriented 
in the axial and circumferential directions respectively and the zjco- ordinate 
directed radially outward (see Fig. la). Thus the axial, circumferential and 
radial deflections at any point on the shell are written as u(x,y,t) , v(x,y,t) 
and w(x,y,t), while the stress field is expressed as an AIRY stress function 
f(x,y,t) such that 


Sx 2 



2.1 



N = ha 

X X 


2.2 


3 



3-4- = -N = - hT 

dxdy xy xy 


2.3 


(h is the wall thickness) 

The use of the AIRY stress function here serves to satisfy the two- 
dimensional equilibrium conditions, 


dN dN 

ox dy 


2.4 


dN dN 

y + xy _ n 

5T 5T'° 

For the case of a thin-walled shell with shallow imperfections, that 
is, a shell satisfying the conditions: 

h « R 

| w(x,y) | « R 


dw 


dw 


si 

i « 1 1 

3p 1 

C< 1 


i 2 v 


^2 
O w 


E 1 57 

or 

a? 


< 0 ( 1 ) 


2.5 


then the following equations hold; 
Compatibility Equation: 


1 4 1 2 2 

V f - — (w, - w, ) -w , + w , +w, w, - w, -w, =0 

Eh R ’XX — XX xy — ’xv ’re ’vv — ’yy — ’w 


xy yy -’xx -’yy 


2„6 


Equilibrium Equation: 
,3 


Eh~ 


12(l-v^) 


V (w-w) + i f -w, f, -w, f, + 2w, f , + phw , = 0 

— R xx ’xx ’yy ’yy ’xx ’xy ’xy tt 


2.7 


where 


’ XX dx 2 


etc . 


Ji a 4 
v = ^ + 


2d 4 , d 4 


dx n dx 2 dy 2 dy 4 
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2.2 Deflection Modes 


The radial deflection function assumed in this analysis is of the form: 


w(x,y,t) = W Q (t) + W^(t)cosKx + W 2 Q(t)cos2Mx + ^ ( t ) c os Mxco sNy + WQgCOs2Ny 

2.8 

while the imperfection function has a similar form: 

w(x,y) = W^cosKx + + W^cosMxcosNy + W^gCos2Ny 2.9 


The first two trigonometric terms of these functions both represent 
axisymmetric radial wall displacements. While these terms are identical in form, 
the provision of two independent wavelength parameters is important in that, 
when investigating the behaviour of shells with : sinusoidal axisymmetric imper- 
fections, it is possible to consider either the case where the axial wavelength 
of the w cosMxcosNy deflection mode is directly related to the wavelength of 
the imperfection ■ (that is, where the imperfection is expressed as W 2 qCos2Mx) 

or, where the wavelength is independent (where the imperfection is described 
by W^cosKx) . The advantage of this flexibility will become more apparent later 
whenthe case of shells with imperfections whose wavelengths differ greatly from 
the critical wavelengths is considered. 

Through the choice of appropriate values for the circumferential wave 
number "N" it is possible to provide geometric continuity with respect to 
radial deflections at y = 0, 277R. However, the deflection function used here 
will not satisfy experimentally realistic end boundary conditions at x =o° 3 L. 
Nevertheless, it has been shown both experimentally 12 and theoretically that 
the effects of clamped end constraints are generally localized to the regions 
adjacent to the ends and that for shells with reasonable length to radius ratios 
(greater than about 2) the effect on ultimate buckling loads is relatively 
small. Accordingly, for the present theoretical study, the effects of end 
constraints are ignored and an "infinitely long" model is considered. 

2.3 Stress Function 


If the assumed displacement function is substituted into the compati- 
bility equation, the equation is identically satisfied by an AIRY stress function 
of the form 

f(x,y,t) = F Q (t) + F^(t)cosKx + F^( t) cosKxcos2Ny + F 2 (t)cos(K+M)xcosNY 
+ F^(t)cos(K-M)xcosNy + F 2 Q(t)cos2Mx + F^^( t)cosMxcosNy 
+ FQ 2 (t)cos2Ny + F 22 (t)cos2Mxcos2Ny + F^^(t )cos3MxcosNy 
+ F^(t)cosMxcos3Ny 2.10 


5 



where 


F = -i £ hy 2 -k £ hx 2 
o 2 x J 2 y 


F l = 


(w -w ) 

Eh v 1 -V 
R 


K 


F = 


F = 
3 


F 4 = 


_ Fh 4F — fw W - W W ) 

( K 2 +4N 2)2 ( 1 02 -1 K)2^ 

K V 

- Eh 2 — "2 ( W L W !r W L W,,) 

2 ( ( K+M) 2 + N 2 ) 2 1 11 1 11 

K 2 / 

— - g — 2“2 (w i w ir \ Hii } 

2((K-M) + it) 1 11 1 LL 


Eh 


F = 
*20 


11 


Eh 


[" (W -W ) + ■ (w 2 -W 2 ) 

1 - ' on Ion' o ' 11 — FI' 


16M 4 1 R 20 - 20 

M 2 


— — (W -W )+ 2M 2 N 2 (W 

( M 2 +W 2 ) L R 11 “ U 


w - w w 
20 11 -20 11 


+ w w - w w ) 

11 02 -11 02 ; 


02 


22 


31 


13 


Eh 

l6N c 


(W 2 -W 2 )' 

2 K ll - IV 


2.11 


EhM 2 W 2 ( ^ 20 ^ 02 ~ —20^02^ 
(M 2 + N 2 ') 2 

Eha^K 2 (« 20 « n - w go w u ) 
(9M 2 + N 2 ) 2 
Eh^N 2 (W U W 02 - 
(M 2 + 9N 2 ) 2 


Here £ x and £ refer to uniform compressive stresses. A full description of the 
derivation of these terms is contained in Appendix A. 

2.4 Mode Equilibrium Equations 

The five-mode displacement and its associated stress function may now 
be substituted into the equilibrium equation. The displacement function chosen 
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here represents only an approximate solution to the equilibrium and compati- 
bility equation set so that the equilibrium equation will not be completely 
satisfied. The error involved may be minimized, however, by satisfying the 
set of equations produced by applying a Galerkin procedure to the left side 
of the equilibrium equation. 

Thus, if E[w(x,y,t) , f (x,y,t) ] =0 is the equilibrium equation, 
then E(W ...,F ...) is unequal to zero by some error. This error is minimized 

by satisfying the equations: 



.) dxdy = 0 


.) cosKx dxdy = 0 


2.12 


2.13 


p 27® p L 
J o o 

n 27® p L . 

O J o 

27® p L 


piim pi 

J o J o 


E(W , ... F q ,...). cos2Mx dxdy = Q 


E(W q , . . ,F q , . . . ) cosMxcosNy dxdy = -0 


E(W , . . .F , . . .) cos2Ny dxdy = 0 


2.14 


2.15 


2.16 


A): 


These equations become, when the integrals are evaluated (see Appendix 


phW = 


. Z h 

phW o = 

? 2 2 

Fh 3 U KF 1 KN r 

K (W -W )+ i + 

2s K 1 -1' R 


2.17 


12( 1-v ) 


T 


VY + 8w 02 f 2 


+ K W Z h 
1 x 


2.18 


P hw 20 . = 


4 Eh 3 


3(l-v £ ) 


14* (W -W ) + ^-.F + 

2s V 20 —20' R 20 


W 11^ F 11 + F 31^: + 


8w F 

02 22 


+ 4m 2 W. 0 >Z h 
02 i x 


2.19 
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phW 


11 


Eh --~ N 2 ) 2 (W -W ) + ^ F 

12(l-v 2 ) 11 11 R U 


2 M 2 N 2 

K 2 ^ 


W (F + F )+W (F +F )+W (F + F ) 
2(T 11 31 ; IF 20 02 J 02 V 11 13 ; 


W, (F + F.) + M 2 ^. Z h + N 2 W 1 Z h 
l 3 ^ 11 x 11 y 


2.20 


phW 02 = - < W 02-V + ^ 


W (F + F ) + 8 w F 

ir 11 13 ; 20 22 


+ 2 K 2 N 2 W l F 2 + 4 n 2 W Q2 Z y h 


2.21 


The 'dot J notation on the left side of these equations refers to differentiation 
with respect to time. 


The set of equations 2.17, 2.l8, 2.19, 2.20 and 2.21 govern the radial 
deflection of a cylindrical shell wall subjected to time-dependent stresses. The 
static buckling equations are essentially identical to equations 2.17 to 2.21 
except, of course, that the time differentials would be equal to zero. The static 
equations are simply a set of non-linear simultaneous algebraic equations instead 
of a group of coupled, time-dependent differential equations. Various methods have 
been devised to handle the general static equations. These generally involve either 
numerical solution schemes or simplifying assumptions. Although a full solution 
of the static equations is beyond the scope of this discussion, it is quite instruc- 
tive to consider the very simple case of the buckling of a perfect shell in a 
single mode. Thus, if it is assumed that all time derivatives, all imperfections 
and each of W^(t), ^^(t), W^^t) and Z are set to zero, then Equation 2.20 may 
be re-written as: ^ 


Eh~ 


12( 1-v ) 


(J^+ N 2 ) 2 


W + 
11 


M 


(m 2 +n 2 ) 2 


Eh TT 

S w n 

R 


k h 

+ (M + IQ 

IS 


EhW 3 = A n Z h 
11 11 x 


2.22 


If we further restrict our attention to the case of small, finite deflections (that 
is, if we consider the situation at the point of buckling), the equation becomes 


Z = E 
x 


L 12(l-v 2 ) 


(M^ N 2 ) 2 


M 2 


M 


(m 2 +n 2 ) 2 r 2 -J 


2.23 


This is an expression for the theoretical buckling stress of a perfect 
shell loaded axially. It is quite simple to show that the minimum value for this 
stress occurs when the wavelength parameters obey the relationship 


(M 2 


+ N 2 ) 2 _ [12( 1-v 2 ) ] L / 2 

2 hR 


M 


2.24 
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This minimum stress, written as, 


. Eh 

R [ 3(l-v 2 )] 1 / 2 


2.25 


is termed the classical buckling stress. 

Although Equation 2.24 does not imply any specific values for either 
M or N, it is worthwhile considering two special cases, namely the case where 
M = N so that: 


M = N 


3(l-v 2 ) 
L 4 h 2 r 2 


1/4 


2.26 


and the case where N = 0 so that: 



2 - V 4 

( 1 -v 2 ) 

2 2 
h R 


2.27 


It may be seen from Equations 2.26 and 2.27 that the wave number of the optimum 
axisymmetric mode is just twice the axial wave number of the optimum "checker- 
board" type asymmetric mode. The two optimum modes may thus be written as 
W2 o cos 2Mx and W^cosMxcosNy with M taking the value defined in Equation 2.26. 

This value is referred to as the critical axial wave number. 


2.5, The Closure Condition 


In addition to providing continuity in w, u and v at y = 0, 277R, it 
is necessary to satisfy the closure condition, that is, to assure that the cir- 
cumferential distance around the shell at any point on the axis is just sufficient 
to close the shell without any overlapping or gaps. This condition may be ex- 
pressed mathematically as 



where 



cr 

_y 

E 


vcr 

x 

E 



w-w 

R 


2.28 


2.29 


The first two terms in this expression represent elastic strains due to 
local stresses in the material, the third term accounts for local slopes or 
folding effects and the last term takes into account the change in circumference 
arising from changes in the radius of the cylinder. 


Recalling Equations 2.1 and 


CT Y (x,y) 

A 


1 

h 


d 2 f(x,y) 

dy 2 


9 


5 



tfyU.y) = 5 


= i d f(x,y) 
" h Sx 2 


w(x,y) = W q + W^cosKx + W 2 qCos2Mx + W^cosMxcosNy + W^cosNy 
w(x,y) = W^cosKx + W^^cos2MIx + W^cosMxcosNy + WQ 2 cos2Ny 


the closure condition may be expanded and eventually reduced to the expression 


O O p p p 

Z - vZ + E 7 T- (w - W ) + N (W - W ) + 
y X [ f ' 11 ; lr ^ 02 - 02' 


W -i 

o 

R J 


In the absence of any external forces such as external or internal 
pressure, the uniform circumferential stress £ is directly related to the uni- 
form radial inertia by Equation 2.17 so that it is possible to write: 


v£ h 

PhW = — f- 
o R 


Eh N' 

T F 


(W 2 U - W 2 L1 ) + N 2 (W : 


2 - W 2 ) + 

02 - 02' R 


2.6 Axial Compression 


The axial compressive deflection of the shell ends may be written as 


or, on the average, 


6 (y) = - f u > x ( x >y) 

n 


- p 27 fR n L 

5 x = - m I I u ’x (x ’ y ) 
'j o ^ o 


Here, as in the case of circumferential extension, the axial strain term 
u , (x,y) is a function of both elastic stresses and radial displacements, 

X 

that is : 


(x M Zl w ! r 

’x^ X ’ y ' E &'>■ ; E ' “ 2 L \d x 


2 / dw \ 2 


With some algebraic manipulation then: 
6 £ v£ ? 


r = r - ¥ + r + r (w \r^u> 


Equations 2.l8 to 2.21 and 2.31 comprise a coupled system of differential 
equations with time-dependent solutions W (t), W^(t) ,W 2Q (t) jW^^t) and W Q2 (t) 
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while Equation 2. 35 relates the load parameters L (t) and £ (t) to the axial 
compression. Together, these represent the governing equations of a cylindrical 
shell with regular sinusoidal imperfections subjected to a time-dependent axial 
compressive load. 

2„7 The Four Mode Solution: A Special Case 

As mentioned in Section 2.2 above, the five-mode deflection function 
used in this analysis provides an important degree of flexibility in that it 
allows for complete independence between the wavelengths of the axisymmetric 
imperfection mode W^cosKx and the other deflection modes. Nevertheless, a care- 
ful examination of Equations 2.l8 to 2.21 (with special attention to the stress 
function coefficients) will reveal that the greatest degree of coupling between 
the various modes occurs when K = 2M„ This is equivalent to saying that the two 
axisymmetric modes W^cosKx and W^qCos2Mx have merged and that the deflection 

function has been reduced to the form 


w(x,y,t) = W Q (t) + W2Q(t)cos2Mx + W^(t)cosMxcosNy + WQg(t)cos2Ny 2.36 

The governing equations for this particular case are identical to those derived 
in this chapter except that W^(t) and are set to zero throughout. 

2.8 Solution 


The solutions for the deflection mode equations were obtained by 
numerically integrating Equations 2.l8 to 2.21 and 2.31 on, an IBM 1130 digital 
computer using a fourth-order Runge-Kutta algorithm. The computer program used 
here was written to describe the behaviour of a shell compressed axially at a 
pre-determined rate. Values of the radial deflection coefficients (W (t), 
W^(t), W 20^ 3 e ^ c *)> the axial load and non-dimensional axial load* were 

computed and printed for successive increments of "time" as the shell was 
compressed. A listing of the computer programs and several sample print-outs 
are contained in Appendix B. 

The computer runs present what is essentially a time history of 
shell behaviour under axial compression. The buckling points were normally 
identified by a sharp drop in axial load accompanied by a rapid increase in 
magnitude of the radial wall deflection coefficients. Figure 2 shows the 
results of a typical run in graphic form. 

When this mathematical model was used to simulate the buckling be- 
haviour of a real shell specimen, the parameters h, R and L were assigned 
values corresponding to the thickness, radius and length of the shell, 
respectively. Likewise, the axisymmetric imperfection amplitude and wave- 
number parameters and K (or W^q and M) were matched to the geometric 

characteristics of the real shell imperfections. For the theoretical repre- 
sentation of both "perfect" and imperfect shells, a very small initial imper- 
fection in the W^°cosMxcosNy mode was specified. This term was not meant as 

* The non-dimensional load was obtained by dividing the actual load by 
the classical static buckling load of the perfect shell with equivalent 
dimensions . 
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a direct representation of a specific physical feature of the shell but rather 
as an initial perturbation which would trigger the growth of the radial deflec- 
tion modes o In the absence of such a perturbation, the theoretical model, being 
mathematically perfect and immune to any disturbances, would be unable to respond 
to instabilities,, 

In contrast with the parameters such as h, R, L and K, which are fixed 
by the geometry of the shell configuration, the buckling mode wave numbers M and 
N are not based directly on shell geometry but rather chosen in such a way as to 
minimize the buckling loads. For perfect shells and shells with imperfections 
with wavelengths close to the critical value, it was generally found that the 
optimum wave numbers could be estimated quite closely by assuming that M = K/2 
and calculating the value of N from Equation 2.24 or 2.26. The exact value of 
the optimum wave numbers was found by varying these estimates slightly and re- 
computing the buckling loads until a true lower bound was established. 

It should be noted that all solutions were obtained for an infinite 
shell model, neglecting the effects of boundary conditions. Based on available 
work for static buckling, the load reductions associated with a clamped end con- 
straint (which was the case for the test models) are of the order of 10$. 
Furthermore, the equilibrium equations did not include in-plane inertia effects. 
However, up to the buckling load, the prebuckling deformation is axisymmetric 
and thus the circumf erential inertia term is zero. On the other hand, the axial 
inertia component is composed of contributions from the rate of the change of 
the closing velocity (V) and prebuckling radial acceleration terms. After an 
initial loading phase, the closing velocity is approximately constant and thus 
the axial inertia can be considered small relative to the radial inertia. 

3. EXPERIMENTAL TECHNIQUE 


As mentioned in the introductory chapter, a basic objective of this 
study has been to produce reliable experimental data on the behaviour of near- 
perfect and imperfect cylindrical shells subjected to dynamic loading. This 
implies three important steps, namely, a method of manufacturing shell specimens 
to very precise tolerances, a means of accurately measuring the shells and 
characterizing their imperfections (including both the imperfections machined 
into the shells for experimental purposes and those arising inadvertently from 
the manufacturing technique) and a method of testing the shells under controlled 
and reproducible conditions. In this study, several modifications of the spin- 
casting technique were used to produce shells with a variety of pre-deter mined 
imperfection geometries. Each shell was carefully measured and then tested on 
a special dynamic loading machine designed by the author specifically for these 
experiments (see Figure 6). 

3.1 Manufacture of Shell Specimens 


The material from which the experimental shell specimens were manu- 
factured was an epoxy plastic mixed from Hysol type XC9~C4l9 resin and type 
356l hardener. This particular material was chosen for its consistent mechani- 
cal properties and because it could be mixed, cast and cured at normal room 
temperatures (thus simplifying the manufacturing procedure) . Prior to the 
manufacture of the shells, an extensive series of tests were carried out on 
samples of the epoxy material to determine the stress-strain characteristics 
under static and dynamic loading. These experiments are discussed in Appendix 
C, 
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The shell specimens were prepared using a spin-casting or centrifugal 
casting technique. In its most basic form, this technique involved pouring a 
small quantity of liquid epoxy , into a cylindrical mold which was spun at high 
speeds (approximately 1750 rpm) so that the material was held in a thin layer 
on the inner surface of the cylinder as it set (see Figure 3)« By coating the 
inner surface of the mold with a liquid releasing agent prior to casting the shell 
layer, adhesion between the epoxy layer and the mold surface was eliminated and 
it was possible to slide the hardened shell out of the mold intact. With precise 
balancing and centering of the mold, this method is capable of producing extremely 
regular shells of essentially constant wall thickness. For example, for the two 
near-perfect shells used in the experimental part of this work, the maximum imper- 
fections and thickness variations in the shell walls were less than 5% of the 
mean wall thickness. Geometric data on these two shells is contained in Table la. 

Two series of shells were produced with major axisymmetric geometric 
imperfections. For the first series, the shells were produced with regular 
sinusoidal profiles machined into the inner surface of an otherwise regular and 
uniform shell. This gave the effect of a regular axisymmetric imperfection of 
the median surface of the wall combined with a sinusoidal variation in wall 
thickness (see Figure lb). The manufacture of these shells began with. the casting 
of a near-perfect shell of slightly greater thickness than required for the final 
imperfect shell. Before this shell blank was removed from the mold, however, the 
mold and shell were mounted on a lathe and carefully centered and circularized. 

The inner surface of the shell was then machined to the desired profile by means 
of a tool-bit controlled by an hydraulic tracer device. This apparatus reproduced 
the profile of an external template onto the inner surface of the shell. The 
development of this system and the preparation of the templates is described in 
some detail in Reference 12. 

The final series of shells produced for this study had matching axis- 
ymmetric profiles on both the inner and outer surfaces, giving effectively a 
constant thickness shell with axisymmetric imperfections of the wall mid- surface. 
These shells were produced using a discardable liner technique. In this method, 
the mold was prepared as described above. However, before the shell was cast, 
releasing agent was applied to the mold and a relatively thin initial layer of 
epoxy was cast in place. When this layer had set, the desired profile was cut 
into its inner surface, exactly as had been done with the previous inperfect 
shells and the profiled surface was coated with releasing agent. Next, the shell 
layer was cast and the inner surface of this layer was subsequently machined to 
the same profile. The shell and the profiled liner were then removed together 
and the liner was split away from the shell (see Figure 4). A total of five shell 
specimens was produced in this fashion, four with regular sinusoidal imperfections 
of various amplitudes and wavelengths and a fifth in which the imperfection profile 
consisted’ of a combination of three superimposed sinusoids. This last shell, 
because of the irregularity of its imperfection profile and because each sinusoid 
contained in the profile existed only as a component of a more general wave-form , 
was considered to have a quasi-random imperfection profile. 

The ends of each shell were cast into fitted grooves on aluminum end 
plates using the same epoxy plastic from which the shells had been manufactured,. 
This procedure provided a very firm and rigid mount. 

3.2 Measurement of Shell Geometry 


A special rig was developed to facilitate the accurate measurement of 
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the geometry of the shell specimens (see Figure 5) = The basic measuring elements 
in this rig were a pair of variable transformer -type linear displacement transducers 
used as contact probes. These probes were mounted on an open, movable frame so 
that they simultaneously contacted the shell wall at the same point on its inner 
and outer surfaces. A motorized lead-screw arrangement moved the probes along 
the shell wall, parallel to the shell axis, at a steady rate. The contact probes 
were connected to demodulators which provided a D.C. voltage signal proportional 
to the displacement of the probe heads. These signals could either be monitored 
directly (to provide profile measurements of the inner and outer surfaces) or 
added electronically by means of an operational amplifier circuit to provide a 
differential signal proportional to the wall thickness. 

The average shell wall thickness for each shell specimen was measured by 
integrating the thickness outputs for a number of axial traverses of the measuring 
rig along the shell. These integrals were converted into thickness deviation 
values which were used to correct the thickness readings made close to the ends 
of the shell with a normal screw- type micrometers. 

For the shells with imperfections cut into the inner surface only, 
the median surface imperfection profile was related to the thickness variation 
of the shell so that the imperfection amplitude could be measured quite easily 
from an X-Y plotter trace of the thickness profile. In the case of shells with 
matching profiles on both surfaces of the wall, it was necessary to measure one 
of the wall profiles independently as well as measuring thickness. For the 
thickness measurements, the arrangement of the two probes and the analogue adder 
tends to cancel out any effects arising from minor misalignment of the shell with 
respect to the measuring rig or any distortion of the shell under probe forces. 

When monitoring the output from a single probe, however, it was absolutely essential 
to have the shell and measuring rig precisely aligned and adjusted for minimum 
contact pressure. The traces recorded from the shells used in these experiments 
are shown in Figures l4 to l8. 

3.3 Dynamic Testing Machine 


A special pneumatically driven testing machine was designed and built 
for the dynamic shell tests (see Figure 6). 

The basic objective in designing this machine was to provide well- 
controlled dynamic loads at a variety of loading rates and to accurately measure 
the applied load and axial end-deflection of the shell during the loading-buckling 
cycle. Figure 7 shows a schematic drawing of the load cylinder of the testing 
machine. In operation, both the upper and lower chambers of this cylinder were 
filled with compressed gas normally to a pressure between 50 and 100 psi. If the 
pressure in the upper chamber were decreased by allowing the gas to escape, the 
resulting pressure differential across the piston caused a compressive force on 
the shell. The rate of the pressure drop in the upper chamber was controlled 
by the size of the venting orifice and by the magnitude of the original pressure. 
Thus, quasi- static loading rates were obtained by relatively small openings of 
the vent valves in the supply control plumbing, while relatively high rates 
were obtained by allowing the diaphragm in the main vent to rupture. A number 
of chokes were manufactured with different size orifices to provide variations 
in the venting (and hence loading) rates. The material used for the diaphragm 
was thin (.002 inches) mylar sheet. 

Due to the high pressures involved, the flow in the venting orifices 
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was supersonic over most of the venting cycle so that the pressure in the upper 
chamber would theoretically vary with time according to 

P(t) = P o (l - e‘ kt ) 3.1 

where P q is the initial pressure and k is a constant depending on the orifice 
size, chamber volume and gas properties. 

By selecting initial pressures sufficiently higher than the pressures 
required to buckle the shells, it was possible to obtain pre-buckling loading 
rates that were approximately constant with time. 

For all tests, a mechanical constraint system was provided to limit 
the axial compression of the shell. Thus, while the shells were allowed to 
buckle completely, large irreversible displacements were avoided and consequently 
the shells could be tested many times and yet continue to yield highly reprodu- 
cible results. 

3.4 Instrumentation 


The testing machine was fitted with instruments designed to measure 
and record the axial loading and deflection of the shell during the loading- 
buckling cycle (see Figure 7)° 

The axial displacement was measured with an MTI model K 45A "FOTONIC 
SENSOR". This instrument consisted essentially of a light source and photo-cell 
masked from each other and mounted adjacent to the end of a divided fibre-optic 
light conducting probe. The opposite end of this probe was mounted close to the 
underside of the lower platen of the testing machine. In operation, light from 
the source was piped down the fibres of the probe and directed onto a small 
piece of reflective aluminum foil fastened to the platen. Thus, as the platen 
was moved with respect to the probe, the amount of light reflected back on 
the end of the probe and transmitted back to the photo-cell changed. Through 
the use of a suitable calibration curve, the output of the photo-cell could be 
related to the displacement of the moving platen. 

A Kistler type 902A load washer was mounted between the piston shaft 
and the lower platen of the testing machine to measure the axial forces exerted 
on the shell. For dynamic loading it was necessary to correct the forces measured 
at this point for the inertia forces arising from the acceleration of the_platen 
and shell end plate. For this purpose, a Kistler type 801 accelerometer was 
mounted at the centre of the lower platen*. The displacement and net force sig- 
nal (consisting of the load washer output less the accelerometer output, suitably 
attenuated) were displayed on a dual-beam oscilloscope (Tektronix type 565) and 
the traces from each run were recorded photographically. Typical oscilloscope 
photographs are shown in Figures 8 and 9« The experimental buckling loads were 
measured directly from these photographs, with a small adjustment being made to 
account for the static pre-buckling load due to the weight of the upper end 
plate. The platen closing velocity "V" was determined from .the slope of the 
displacement -time curve immediately preceding buckling. 


* The accelerometer, which consisted basically of a sprung mass, exhibited some 
high-frequency ringing which introduced unwanted noise into the final force 
signal. It was found necessary to provide a low-pass electronic filter (minus 
3 db at 5000 Hz) between the accelerometer and the oscilloscope used to monitor 
the signals . 
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In addition to the normal tests, several runs were made in which various 
shells were photographed with a high-speed movie camera during loading and buckling . 
The epoxy plastic material from which the shells were manufactured was translucent 
and birefringent so that rather than photographing the shells directly, it was 
convenient to use photo- elastic techniques and to photograph the 4 5 isoclinic 
patterns as they arose in the shell. This was achieved by illuminating ;the shell 
model with a powerful quartz-halogen flood lamp, the light from which passed 
through a polarizing filter mounted with its polarizing axis inclined at 45 to 
the axis of the shell (and hence at 45 to the principal stress direction) . The 
light was reflected from a thin layer of fine aluminum powder deposited on the 
inner surface of the shell and thence passed through a second polarizing filter 
mounted' on its polarizing axis perpendicular to the axis, and onto the camera 
lens. The camera and the flood lamp were located as close together as possible 
and on the same level as the centre of the shell so that the direction of the 
incident light and the optical axis of the lens were as close as possible to 
being perpendicular to shell wall. The camera used was a HYCAM l6 mm. model 
with a 4 00 foot film capacity. This instrument was capable of operating at a 
rate of up to 5000 full frames per second. Because of the short exposure times 
resulting from the very high framing rate and the light losses due to the 
polarizing filters, it was necessary, to use the maximum possible illumination 
intensity and extremely sensitive film materials. Due to the large amount of 
heat associated with a high intensity light beam, it was necessary to use several 
heat absorbing (infra-red) filters immediately adjacent to the lanp to avoid 
damage to the first polarizing filter and to minimuze heating of the shell 
specimen. Further, a variable transformer was used in the power supply to the 
lamp so that the amount of illumination could be regulated, full intensity being 
used only during the actual filming. Eastman Kodak 2475 and 2485 high-speed 
recording films were used in these experiments. These films were processed in 
highly energetic developers to achieve full contrast. Normally full 400 foot 
spools of film were used to insure that there would be adequate time for the 
camera mechanism to accelerate to its full framing rate. 

Through the use of the photo-elastic technique and high speed photo- 
graphy, ’ it was possible to observe, in detail, the development of internal 
pre-buckling strains in the shell walls. A more complete description of the 
photo-elastic method is found in Reference 29, while a further description of 
the application of the high-speed movie camera in the study of buckling is 
found in Reference 30. 

4. DISCUSSION OF RESULTS 

4 . 1 Preliminary Theoretical Results 


Figure 10a compares a series of axial load versus time curves computed 
for a hypothetical near-perfect shell compressed at different rates. The axial 
loads in this figure are expressed in non-dimensional units obtained by dividing 
the computed loads (P). by the classical static buckling load of the shell 
(P cL ). Similarly, the axial compression rate is given in terms of the non- 
dimensional parameter Vo, equal to the closing velocity of the loading platens 
(v) divided by the acoustic wave speed of the plastic material c = (e/p)V^. 

The same data is re-plotted in a different form in Figure 10b. The curve shown 
here is an envelope of the maximum axial loads and represents the relationship 
between the dynamic buckling loads and the rate parameter Vo. The results indi- 
cate that there is a substantial increase in the magnitude of the buckling loads 
at the higher compression rates. This is due to the effect of inertia, which 



tends to retard the growth of the wall deflections. 

The influence of regular axisymmetric imperfections on dynamic buckling 
loads is shewn in Figures 11 and 12. Figure 11 contains a family of curves 
showing the variation of buckling strength with imperfection amplitude (for 
critical wavelength imperfections) while Figure 12 shows the effect of different 
wavelengths with a fixed amplitude. The results given in these figures indi- 
cate that the effect of imperfections on dynamic buckling loads is qualitatively 
quite similar to their effect on static loads. In both cases, the buckling 
strength may be seriously impaired by the presence of relatively shallow .imper- 
fections. Further, the theoretical results indicate that the effect of sinu- 
soidal imperfections is much more serious when the wavelength of the imperfec- 
tion profile is close to the critical wavelength for the shell in question. 

When computing the theoretical buckling loads used in constructing the 
graphs shown here, the "independent" wave-numbers M and N were in each case 
chosen to provide minimum buckling loads. In the cases where the wave-number 
of the axisymmetric imperfection was taken t*b be nearly equal to the critical 
wave-number, it was found that the four-mode solution (that is, the case where 
K = 2M) provided the lowest theoretical buckling points.. However, for the 
special case where the imperfection wave-number was substantially smaller than 
the critical (that is, the imperfection wavelengths were longer than critical) 
it was found at the higher deflection rates that the alternative five-mode 
solution would yield lower results. This suggests that where the imperfection 
wavelengths are relatively close to the critical, the geometry of the buckling 
deflections will be largely predetermined by the imperfection mode shape, but 
that in cases of highly non-optimum imperfections, buckling mode shapes will 
be essentially independent of the initial imperfections, and unless these im- 
perfections are quite large, they will have a relatively small influence on 
the buckling strength of the shell. 

A basic objective of this study has been to provide a direct com- 
parison of experimental and theoretical results. The dynamic loading configu- 
ration chosen as being most readily adaptable to both experimental realization 
and theoretical modelling was compression at a constant rate of axial deflec- 
tion. This corresponds to an axial deflection function of the form 

6(t) = V t 

Although the dynamic testing machine used in these experiments proved capable 
of providing reasonable constant deflection rates under most conditions, it was 
apparent from the experimental results that at the higher rates, the time re- 
quired to accelerate the moving piston, platen and end plate became significant 
when compared to the total time of the test and that the axial deflection rate 
was considerably lower at the beginning of the test than at the buckling point. 

In order to determine the theoretical consequences of this situation, a number 
of load- time histories were computed using a deflection function of the form 

6(t) = V (l-e" 7t ) t 

The results, shown in Figure 13, indicated that variations in the 
loading rate during the initial stages of the test have very little influence 
on the ultimate buckling strength of the shell. It was concluded, therefore, 
that the dynamic load history of the experimental shells could be adequately 
characterized by a single parameter, that is, the deflection rate measured in 
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the linear region of the deflection-time curve immediately preceding the buckling 
point. 

4 . 2 Measurements of Shell Geometry 


Ten* shells were manufactured for this project. These included two 
near-perfect shells, three shells with axisymmetric sinusoidal imperfections of 
various wavelengths and amplitudes machined on the inner surface only, four shells 
with matching inner and outer surfaces (constant thickness imperfect shells) and 
one shell with matching inner and outer profiles which consisted of a superposition 
of several sinusoids of different frequencies. The characteristics of each of 
these shells are summarized in Table 1. 

Traces of the profiles of the imperfect shells are shown in Figure l4 
to 18. For shells 3 to 5, the traces were obtained from thickness measurements 
and represent the median surface profiles. Due to the averaging action of the 
probes, any overall distortion of the shell surface away from a perfect cylinder 
defined by the shell ends is ignored. For shells 6 to 10, the outer surface pro- 
file and thickness were measured separately along four equally- spaced generators 
on each shell. In Figures l6 to l8 the upper four traces show the outside surface 
profiles while the lower traces show thickness variations. These traces show that 
the periodic thickness variations of the shells were generally small when compared 
to either the average thickness or to the imperfection amplitude. The thickness 
variations that do exist are primarily due to a slight axial misalignment of the 
inner and outer profiles (that is, a phase mismatch) rather than a difference 
between the amplitudes of the profile amplitudes. The outer surface profiles 
shown in these figures indicate a slight inward bow of the shell walls between 
the end plates. This variation in radius is more apparent than real, however, as 
it arises largely from deflections in the wall under pressure from the outer 
contact probe. 

The profile measurements described here demonstrate that the shell 
specimens used in this study were very uniform in their overall geometry and 
that the random imperfections and thickness variations arising from machining 
errors and other manufacturing faults were generally much smaller than the pre- 
scribed imperfections. 

4.3 Comparison of Theory and Experiment; Dynamic Buckling Loads 

An important feature of the epoxy mixture used in the manufacture of 
the shell specimens is that this material is capable of sustaining quite large 
strains without permanent deformation. Consequently, it was possible to buckle 
the shells repeatedly under either static or dynamic loads without damaging the 
material or altering the buckling strength of the shell. This meant that an 
extensive series of tests at different loading rates could be made with a single 
specimen. The experiments described in this report represent more than two hundred 
buckling tests. This includes 124 carefully controlled buckling strength tests 
(which provided the data points shown in Figure 19 through 27), four tests which 
were filmed with the high speed movie camera, ten tests made with restricted 
axial deflections (see Section 4.5) and nearly one hundred pretest runs made 
during the development and calibration of the testing rig. 

Figure 19 shows the experimental results obtained with two near-perfect 
shells, plotted together with the theoretical buckling load-deflection rate curve. 
These 15 experimental points lie consistently between 5 °lo and 15^ below predicted 
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buckling loads for the entire range of axial deflection rates.. Except for this 
small offset, which may readily he attributed to the effects of end constraints 
and to small residual imperfections, the agreement between theory and experiment 
is very good. 

The results obtained from b5 tests with the three shells with axisym- 
metric imperfections on the inner surface only are shown in Figures 20, 21 and 22. 
The theoretical buckling loads of these shells were computed using the four-mode 
model with initial imperfections in the cos2Mx and W cosMxcosNy modes. In 

each case, the axisymmetric imperfection amplitude (W^q) and wavelength were set 

equal to the amplitude and wavelength of the mid-surface imperfection, while the 
W ^ cosMxcosNy term served as a small asymmetric perturbation. 

When comparing the theoretical and experimental results, it is apparent 
that although the experimental points lie within about 20$ of the theoretical curve, 
these points indicate a more rapid increase in dynamic buckling load with deflec- 
tion rate than is predicted theoretically. In several instances, the experimental 
buckling loads are even higher than those predicted by the somewhat idealized 
theoretical model. An earlier analysis of imperfect shells of this configuration^ 
had shown that for the case of static buckling, the presence of small-thickness 
variations has a relatively minor effect on shell stiffness and stability. For 
the dynamic case, however, the periodic mass distribution arising from the thick- 
ness variations evidently gives rise to inertia effects which significantly increase 
the dynamic strength of the shells. This effect is not considered in the mathe- 
matical model which assumed constant thickness shells. 

In order to eliminate this discrepancy between theory and experiment, a 
further series of shells was produced with matching sinusoidal imperfection pro- 
files on both the inner and outer surfaces so that the walls were of essentially 
constant thickness. The theoretical and experimental results obtained with these 
shells are shown in Figures 23 to 27. A total of 55 data points are reported. 

For the first two shells of this series (shells number 6 and 7 ), the 
wavelengths of the axisymmetric imperfections were quite close to the theoretical 
critical wavelengths. In both cases, the minimum theoretical buckling loads were 
found using the four-mode model and in each case the agreement between theoretical 
and experimental buckling strength is good for the entire range of compression 
rates encountered in the tests. The imperfection wavelength of the third shell 
(shell 8) was approximately half the critical wavelength (that is, the wave-number 
of the axisymmetric imperfection was 'almost twice the critical). Here again, the 
minimum theoretical buckling loads were obtained using the four-mode theoretical 
model and the agreement between predicted buckling strength and experimental 
results was quite good. Shell 9 bad an axisymmetric sinusoidal imperfection 
profile with a wavelength slightly more than twice the critical. Here, the 
theoretical representation of the initial buckling mode was more complicated than 
had been encountered before. It was found that at the lower axial deflection 
rates, minimum theoretical buckling loads were computed using the four-mode 
solution but at the higher rates, lower loads could be obtained when the five- 
mode model was employed. Physically, this implies that the buckling mechanism 
would change when different loading rates were encountered. Thus, under quasi- 
static or gradual loading, the buckling deflection pattern would be dominated 
by the imperfection and would have long axial wavelengths. At higher loading 
rates, however, the optimum buckling deflection modes are relatively independent 
of the pre-buckling deflections (which, because of inertia, are comparatively 
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small) and assume shapes resembling those encountered in the buckling of near- 
perfect shells. The maximum theoretical buckling load at a given loading rate 
would be given by the lower of the two predictions. The theoretical curves, 
together with the experimental results, are shown in Figure 26. 

Except for the two near-perfect shells, the shells used in the tests 
described above all contained regular sinusoidal imperfections. It was recog- 
nized, however, that the imperfections encountered in industrially produced shell 
structures would generally be random in nature with the critical axisymmetric 
wave-form occurring only as a component of a more general profile. In order to 
simulate this situation, a shell was manufactured with an imperfection profile 
that consisted of three superimposed sinusoids. Although this profile was 
basically periodic (see Figure 18) it was similar in nature to the random imper- 
fection case in that the critical wavelength component was contained in a more 
general wave-form. 

The experimental results obtained from 12 tests on this shell with 
quasi-random imperfections are shown in Figure 27. These results are compared 
with two theoretical models. To obtain an upper bound estimate of the dynamic 
buckling strength, theoretical buckling loads were computed for a hypothetical 
shell with a regular sinusoidal imperfection of critical wavelength whose ampli- 
tude was equal in magnitude to the Fourier component of the general profile with 
the equivalent wavelength. As a lower bound, a second hypothetical case was 
considered of a shell with a critical wavelength sinusoidal imperfection having 
an amplitude equal to the root mean square value of the general imperfection 
function. The first of these models represents the case where it is assumed 
that the buckling mechanism is entirely dominated by the critical wavelength 
component and that the other imperfections have essentially no influence on 
the ultimate strength, while the second model represents the opposite extreme, 
that is, assuming that the effect of the various imperfection components is 
cumulative and independent of the actual wave-form. For the imperfection pro- 
file tested, these two models give quite different results and the experimental 
buckling loads all lie between the two bounds. 

14 

It has been suggested that the theoretical buckling strength of 
the axisymmetric root mean squared amplitude model would provide a useful design 
load for structures incorporating thin-walled cylindrical shells with generalized 
random imperfections. While it must be recognized that the results presented 
here are for a single specimen and have no statistical signifance, it has been 
shown that the R.M.S. lower bound does provide a safe but no unduly conservative 
estimate of the dynamic buckling strength of the shell. 

4 .4 Comparison of Theory and Experiment; Buckling Modes 

According to the mathematical model of dynamic buckling developed in 
Section 2 of this report, the collapse of the shells is due to the combined 
influence of the axisymmetric and asymmetric deflections. More specifically, 
the numerical results suggest that for a shell with major axisymmetric imper- 
fections, the radial wall deflections develop in two stages (see Figure 8). 

During the initial or pre-buckling loading, there is a gradual increase in the 
amplitude of the axisymmetric imperfection wave. As the buckling point is 
approached, the asymmetric deflections, which have hitherto been comparatively 
small, begin to grow rapidly so that the buckling mode has major axisymmetric 
and asymmetric components. 
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This buckling mechanism was confirmed from high speed photographs of 
the 45 isoclinics observed through a system of polarizing filters. As mentioned 
in Section 2, this was possible because of the photo-elastic nature of the epoxy 
plastic. A series of frames from a film taken with the high-speed move camera 
during a dynamic test of shell number 3 are shown in Figure 20. This sequence 
includes the final stages of loading, the buckling event and the early post- 
buckling period. The pictures illustrate that during the pre-buckled stages, 
the principal observable deformations are in fact axisymmetric , but that with 
the onset of buckling, a regular 'checker-board' type asymmetric pattern appears. 
Following buckling, the isoclinic patterns begin to shift their shape and 
eventually the diagonally oriented diamond shaped pattern typical of buckled 
shells emerges. Measurements made from the photographs showed that the circum- 
ferential and axial spacing of the asymmetric pattern corresponded fairly well 
(within about 30 ‘jo) to the circumferential wavelengths of the theoretical 
buckling mode. 

It will be recalled that during the discussion of the dynamic buckling 
of shell 9 3 it was hypothesized that at higher deflection rates the buckling 
modes of this shell would be largely independent of the initial imperfection 
wave. A number of frames from a film taken during a high-speed run with this 
shell are shown in Figure 29. These photographs show that while the axisymmetric 
pre-buckling deflections have the same wavelength as the initial imperfection, 
the wavelengths of the buckling deflections are generally shorter and correspond 
quite closely to those which would be expected for a near-perfect shell. 

A quantitative comparison between experimental and theoretical wave- 
forms is given in Table 2. 

4.5 Shell Stability Under Quasi-Stepwise Loading 

The dynamic loading configuration and dynamic buckling loads considered 
so far in this analysis are of primarily academic interest. The designer is 
not particularly concerned with the conditions under which buckling occurs but 
rather with the loads which may be imposed without incurring instability. 

A typical dynamic loading situation of significance to the engineer 
would be the case where an axial load increases rapidly up to a value just below 

the critical and then remains at that level. This is equivalent, for example, 

to the case where a rocket engine fires and quickly achieves a thrust nearly 
equal to the static buckling strength of the fuselage supporting it. It is quite 

clear that if the ultimate axial load on the shell exceeds the static critical 

load, the shell will eventually buckle, regardless of its earlier loading history. 
What is of interest here, however, it whether or not the momentum of the dynamic 
pre-buckling deflections would be sufficient to carry the shell into its buckled 
state even though the maximum load remains below the static stability limit. This 
problem has been discussed a number of times in the literature (see References 
l8, 19, 20 and 23) . Each of these studies has dealt with the case where the 
axial load was applied in a stepwise fashion, that is, where the axial load in- 
creased instantaneously from zero to its maximum value, and each has shown that 
for imperfect shells, the dynamic stability limit lies somewhat below the static 
critical strength of the shell. Due to limitations in the experimental apparatus, 
it was not possible to physically duplicate the stepwise loading case. However, 
the situation was approximately simulated on the dynamic testing machine by 
blocking the travel of the loading platen to restrict the axial deflection (and 
hence axial load) of the shell. At the highest piston velocities available with 
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the testing machine, the rise-time of the axial load was about 1.5 milliseconds. 
Theoretical results computed for this quasi-stepwise deflection-time configuration 
indicated that at experimentally realistic deflection rates, the near-perfect 
shells should remain unbuckled if the maximum axial load does not exceed the 
static critical load. These experimental results were not entirely conclusive. 

In about 5 Qfj 0 of the runs made with the testing machine adjusted so that the axial 
travel was slightly less than that required to buckle the shell statically, the 
shell remained unbuckled. This suggests that the dynamic stability limit is 
roughly the same as the static buckling strength. A more precise experimental 
investigation of this probelm is beyond the capabilities of the present experi- 
mental apparatus. 

5. CONCLUSIONS 


The experimental data presented in this report indicate that the experi- 
mental techniques and shell manufacturing methods adopted in this study provide 
the most reliable and highly reproducible means of determining the dynamic axial 
buckling strength of thin- walled cylindrical shells hitherto available. Further, 
the cldse agreement between experimental and theoretical results for both buck- 
ling strengths and buckling modes shows that the mathematical model chosen in 
this analysis adquately accounts for the principal mechanisms governing the 
dynamic behaviour of these structures. Specifically, it may be concluded that 
the buckling deflection modes of real shells are similar to those employed in 
this model, that the change in buckling stiffness with compression rate is 
principally due to the effect of inertia limiting the rate of growth of radial 
deflections and that the dynamic buckling strength of shells is seriously 
impaired by the presence of geometric imperfections, especially where these 
imperfections are regularly spaced at intervals close to the wavelengths of 
important buckling modes. 

These conclusions have several important implications for the struc- 
tural engineer working with designs involving shell- type members. It is very 
clear that the amplitude and form of the wall imperfections have a major effect 
on the buckling strength of shells under dynamic as well as static loading. The 
changes in slope of buckling load versus loading rate curves indicate that the 
imperfection wavelengths have an especially notable influence on dynamic shell 
response since the wavelengths tend to dictate the deflection modes and hence 
influence the fundamental frequencies of the deflections. 


Although the range of dynamic loading configurations and shell 
geometries considered in this study was restricted by limitations in the 
experimental apparatus, the general success of the mathematical model suggests 
that an analytic approach based on this model could be used very profitably 
to explore a broad spectrum of problems in shell design for dynamic load situa- 
tions . 

O /T 

It has been pointed out by Professor J. Hutchinson 0 that the choice 
of the loading rate parameter, V/c, does not permit the results, as plotted, 
to be compared to other cylinder configurations. The suggested non-dimensional 
parameter employs the period of vibration of the classical "square-wave" mode 
(t(1) = 2W2 R/c) and the time required for the shell to reach the critical 
"static" strain for buckling 
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in the form: 


T (l) 

^27 


2W6(l-.v 2 ) 2 R 2 

hL 


(V/c) 


It would appear that this parameter is more suitable, taking into account both 
the shell geometry and dynamic effects. 
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APPENDIX A: 


I: DERIVATION OF AIRY STRESS FUNCTION 

As explained in Section 2.1, the compatibility equation for an im- 
perfect cylindrical shell may be written as (see Equation 2.6): 
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If it is assumed that the displacement function w takes the form 
w(x,y,t) = W Q (t) + W^(t) cosKx + W^^t) cos2Mx + W^(t) cosMxcosNy 

+ Wyg(t) cos2Ny A. 2 


and that the imperfection function has the form 

w(x,y) = W^cosKx + W^qCos2Mx + W^^cosMxcosNy + W^cos^Ny A. 3 

then the right-hand side of the compatibility equation may be re-written as 

- K 2 (W^-W^) cosKx - 4M 2 (WgQ- cos 2Mx 


follows: 
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When the above equation is fully expanded, there are several higher order triono- 
metric products which may be reduced using the following identities: 

2 2 2 2 

cos Mxcos Ny - sin Mxsin Ny = \ cos2Mx + \ cos2Ny 

cosKxcosMxcosNy = § cos(K+M)x cosNy + \ cos(K-M)x cosNy A. 5 

cos2MxcosMxcosNy = \ cos3MxcosNy + \ cosMxcosNy 
cosMxcosNycos2Ny = \ cosMxcosNy + cosMxcos3Ny 
On expanding and collecting terms, Equation A. 4 becomes: 
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This equation may be satisfied if it is assumed that the stress func- 
tion "f" has the form 

f(x,y,t) = F^(t)cosKx + F 2 (t)cosKxcos2Ny + F^(t)cos(K+M)x cosNy 

+ F^(t)cos(K-M)x cosNy + F 20 (t)cos2Mx + F^t^osMxcosNy 
+ F 02 (t^cos2Wy + F2 I (t)cos3MxcosNy + F^C'OcosMxcosSNy 

Z h Eh 

+ F 22 (t)cos2Mxcos2Ny - y - — ^ x A. 7 

With this premise, the left hand side of the compatibility equation becomes: 


K^F^cosKx + (K 2 + 4N 2 ) 2 F 2 cosKxcos2Ny + ((K+M) 2 + N 2 ) 2 

• F^cos(K+M)x cosNy + ((K-M) 2 + N 2 ) 2 F i+ cos(K-M)x cosNy 

4 .2 2 2 ^ 4 
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+ N 2 ) 2 F^^cos3MxcosNy + (M 2 + 9N 2 ) 2 F^ cosMxcos 3Ny 

+ l6.(M 2 + n 2 ) 2 F 22 cos2Mxcos2Ny A. 8 

For the compatibility equation to remain valid for all x and y, it is 
necessary to equate the coefficients related to each trigonometric expression. 
Thus: 
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cos2Mxcos2Ny coefficients: 
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It is a simple matter to re-write this set of equations as expressions 

for the stress function coefficients F , F^, etc. (see Equation 2.1l). The 

stress function f(x,y,t) is thus completely determined by the displacement 

field w(x,y,t) and uniform stresses (Z and Z ) existing in the shell. 

x y 

II: SOLUTION OF EQUILIBRIUM EQUATION 


The equilibrium equation introduced in Section 2.1 is as follows 
(see Equation 2.7): 
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If the trigonometric series representation of the deflection and stress 
functions given above are substituted into the left hand side of this equation, 
one obtains the expression, 
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If Equation A. 11 is expanded and the trigonometric coefficients collected, 
it then becomes: 
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- 4 ( K+M) 2 N 2 W q 2 F o c o s ( K+M) x cosNycos2Ny-36M 2 N 2 W Q2 F cos3MxcosNycos2Ny 


- 4k 2 n 2 w q2 f 2 cosKxcos 2 2Ny-l6M 2 N 2 WQ 2 F 22 cos2Mxcos 2 2Ny 

- 4M 2 N 2 WQ 2 F^ cos M xcos2 Nycos3Ny - ^y* 1 + phW Q = 0 

R 


A. 12 


Although it would be possible to apply the Galerkin procedure to the 
above expression and evaluate the integrals (Equations 2.12 to 2.15) directly, 
it is possible to achieve identical results with considerably less effort by 
rearranging this expression in a form that avoids multiples of the trigonometric 
terms and by using the orthogonality property: 


27 r 


cos i x cos jxdx=Oifi^j 


A. 13 


This rearrangement of expression A. 12 involves the use of the following trigono- 
metric identities : 


cosA cosB = §oos(A+B)x + |-cos(A-B)x 
sinA siriB = |cos(A-B)x - §-cos(A+B)x 


A.14 


After these transformations have been made and;the terms collected 
according to their trigonometric coefficients, the equilibrium expression 
becomes : 
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k4 (w r V - r F i- K \V - < w n F 3 + w uV 


L 12 ( 1 -v ) 


+ phW_ 


cosKx + 


4 Eh' 


3(l-v 2 ) 


2 

M < W 2CfV - i f- F 20-' ,M2w 20 I: x h - M2l,2( W 11 F 
,3 


+ W 11 F 31 + 8W 02 F 2 2 ) + phW 20 
M 2 


11 11' 


cos2Mx + 


Eh“ 


L 12(l-v ) 


(mV) 2 (w,,- w ,) 


- - F n A ii\ h - < w 20 (f u + f 3! )+ w h (F 20 + V 


+ w 02 <F u + F 13 )) - W V + pl *u 


cosMxcosNy 


4Eh 3 
L 3(l-v 2 ) 


^ (W 02" W 02 ) - i+ N 2 Wno Z xx h - M2 N 2 (W nl (F 11 + F 1Q )+ 8¥ Qn F OQ ) 


02 y 


2K 2 N 2 W 1 F 2 + phW Q2 


cos2Ny + 


ll' 11 13' 

cos3MxcosNy 


20 22 ; 


cos(M+K)xcosNy + 


2 h „ 

-X— + phW =0 
R o 


A. 15 


Applying the Galerkin procedure to this equation is equivalent to simply equating 
the coefficients of the appropriate trigonometric terms to zero. This results 
in Equations 2.17 to 2.21 of t,he main text. 
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APPENDIX B 


COMPUTER LISTINGS AND OUTPUTS 

This section contains a listing of the computer programs used to inte- 
grate the set of time-dependent differential Equations 2.17 to 2.21 and a set 
of sample outputs. The program is written in Fortran IV language and is designed 
to run on the IBM 1130 computer. The subroutine RKGIL feferred to in this pro- 
gram is a fourth-order Runge-Kutta integration algorithm supplied by IBM as part 
of the standard scientific systems package for this machine. Two versions of 
the program and two versions of the subroutine NYST are given here, one for use 
when considering the general five-mode displacement function and one for use with 
the four -mode approximation. 
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FIVE MODE MODEL 


// 

* 

* 

* 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


FOR 

LIST SOURCE PROGRAM 

IOCS (1132 PRINTER. CARD. DISKI 

ONE WORD INTEGERS 

FIVE MODE (THREE WAVENUMBER) RUNGE KUTTA SCHEME 


i-****#*#* 

ITN 

ISHEL 

C.A.B 

L 

R 

H 

Wl.WTO WOO. WOT 

Y ( 1 ) • • • Y ( 5 ) 
Y(6)...Y(10) 

T 

TR 

VO 

DISP 
PAX 
PCR 
SIGEX 
AS . RS . DEL 


************************* 
NO. OF WAVENUMBER COMBINATIONS PROVIDED 
SHELL NUMBER 

AXIAL + CIRCUMFERENTIAL WAVENUMBERS (=K.M.N) 

SHELL LENGTH 
SHELL RADIUS 
SHELL THICKNESS 

IMPERFECTION AMPLITUDES (W1.W20.W11.W02 
MODES RESPECTIVELY) 

RADIAL DEFLECTIONS 

TIME DERIVATIVES OF Y ( 1 ) ETC. 

TIME PARAMETER ( T I ME /ONE-D I M. WAVESPEED) 

REAL TIME 

AXIAL VELOCITY PARAMETER 
(VELOCITY/ONE-DIM. WAVESPEED) 

AXIAL DISPLACEMENT OF SHELL END 
AXIAL LOAD 

CLASSICAL AXIAL BUCKLING LOAD 
AXIAL STRESS 

STARTING POINT. END POINT AND INCREMENT 
FOR RUNGE-KUTTA CYCLE 


EXTERNAL NYST 
DIMENSION Y ( 1 0 ) . P < 5 ) 
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C THREE WAVE-NUMBER RUNGE-KUTTA SCHEME 

READ ( 2 *200 ) ITN 
DO 2 1=1. ITN 
READ (2.201) A.B 

200 FORMAT (110) 

201 FORMAT (2F10.5) 

ISHEL = 9 

C = 2.282 
H=0. 01746 
R=3.900 
Wl=-0. 00294 
WTO=0 « 0 
W00=0. 00001 
WOT=0 .0 
V0=0.0001 
DO 1 J= 1 , 5 
STD=0.0 
P ( 2 ) = VO 
P ( 3 ) = A 
P ( 4 ) = B 

WRITE (3.100) ISHEL. VO. Wl. WTO. WOO. WOT. C. A.B 

100 FORMAT ( 1CH1SHELL N0..I3./.4H V0=»F7.b.4H W1=.F8.5»5H W20=.F8.5.5H 
1 W11=.F8.5«5H W02= . F8 • 5 . 3H K=*F6.3.3H M=.F6.3.3H N= »F6 . 3 . // . 5X »4HT 
2IME.8X.5HWOIT) .7X.5HW1 (T) .6X.6HW20IT) .6X.6HW1KT) .6X.6HW02(T) .7X.4 
3HP ( T ) *5X»8HP(T )/PCL./) 

Y(1)=W1 

Y ( 2 ) = W T 0 
Y ( 3 ) =W00 
Y ( 4 ) =WOT 
Y ( 5 ) =0.0 
Y ( 6 ) =0.0 
Y( 71=0.0 

Y ( 8 ) = 0 . 0 
Y ( 9 ) =0 .0 
Y ( 10 ) =0 . 0 
SIGEX=0,0 
AS=0. 0 
BS=7. 0 
DEL=1 .75 

DO 6 K = 1.1000 

CALL RKGILf Y . 1 0 . AS »BS . DEL . N YST .P) 

T = P ( 1 ) 

S I GEX = P ( 5 ) 

PAX=2477000.*R*H*SIGEX 
PCR=1555000. *H*H 
PND=PAX/PCR 
TR = T /70000 . 

WRITE (3.101) TR.Y( 5) .Y( 1 ) ,Y(2 ) .Y< 3) * Y ( 4 ) .PAX.PND 

101 FORMAT (1X.6IF9.5.3X) .F9.3.3X.F9.5) 

IF (STD-PAX-20. ) 5.1.1 

5 AS=BS 
STD=PAX 

CALL DATSW(l.KK) 

GO TO ( 1.6) .KK 

6 BS=AS+7. 

1 VO=VO+0. 0001000001 

2 CONTINUE 
CALL EXIT 
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SUBROUTINE NYST ( M « T t Y . YDOT .P ) 

SUBROUTINE NYST FIVE MODE MODEL 

F1*F2... STRESS FUNCTION COEFFICIENTS 

NU POISSON'S RATIO 

SIGEX AXIAL STRESS 

SIGWY CIRCUMFERENTIAL STRESS 

DIMENSION YdO).YDOT(lO) .P<5) 

VO=P< 2 ) 

A°P ( 3 ) 

B = P( A ) 

C=2.282 
H»0. 01746 
R=3.900 
Wl=-0. 00294 
WTO=0.0 
W00=0. 00001 
W0T = 0 • 0 
NU=0.40 
Z«2. 52 
L= 1 1 • 0 

F1=-(Y( 1 >-Wl ) / ( R*C*C ) 

F2=-4.*C*C*B*B*( Y ( 1 ) *Y ( 4 ) -W1*W0T ) / ( C*C+4. *B*B ) **2 
F3=- C*C*B*R*IY( 1 ) *Y ( 3 ) “W 1*W00 ) / ( 2 • * ( ( C+A ) **2+B*B ) **2 ) 

F4=-C*C*R*B* ( Y ( 1 >*Y< 3 )-Wl*WOO>/ ( 2.*( <C-A ) »*2+B*B) **2 ) 

F20=- (4./R*( Y( 2 )-WT0) + .5*B*B*( Y( 3 >**2-WOO**2 ) ) / ( 16 .*A*A ) 
Fll=-(A»A/R*lY<3)“WOO>+2.*A*A*B*B*IY<2>*Y<3)-WTO*WOU+Y(3)*Y(4)-WOO 
1*W0T ) ) / (A*A+B*B)**2 
F02=-A*A*( Y ( 3 ) **2-WOO**2 ) / ( 32.*B*B ) 

F31=-2.*A*A*B*B*( Y(2 >*Y( 3>-WT0*W00> / <9.*A*A+B*B>**2 
F13=-2.*A*A*B*B* ( Y ( 3 ) *Y ( 4 ) -WOO* WOT )/(A*A+9«*B*B)**2 
F22=- A*A*B*B* ( Y ( 2 ) *Y ( 4 >-WTO*WOT ) / ( A*A+B*B ) **2 
D I SP= VO*T 

SIGEX=NU*SIGWY+DISP/L-t25*C*C*( Y l 1 ) **2-Wl*W 1 ) -A*A* ( Y ( 2 > **2-WTO**2 + 
10»125*(Y(3) **2 — W00**2 ) ) 

SIGWY=NU*SIGEX-B*B*< Y l 4 > **2-WOT**2+. 125* (Yl 3 ) **2-WOO**2 ) )-Y(5)/R 
YDOT ( 1 >=Y<6) 

YDOT ( 6) =-H*H/(4.*Z >*C**4*(Y< 1 ) -W1 ) +C*C /R*F1+. 25*C*C*B*B* ( Y ( 3 ) * ( F3+ 
1F4 )+8.*Y(4)*F2 )+C*C*Y( 1 )*SIGEX 
YDOT ( 2 ) =Y ( 7 ) 

YDOT ( 7)=-4.*H*H/Z*A#*4*< Y(2 ) -WTO ) +4. *A*A/R*F20+A*A*B*B* ( Y(3 )*(F11+ 
1F31)+8#*Y(4)*F22>+4.*A*A*Y<2)*SIGEX 
YDOT ( 3 ) =Y ( 8 ) 

YDOTf 8 >=-H*H/<4.*Z)*(A*A+B*B)**2*<Y( 3 > -WOO ) +A*A/R*F 1 1 +2 . *A*A*B*B* < 
1Y(2)*(F11+F31)+Y(3)*(F20+F02)+Y(4)*(F11+F13) > +. 5*C*C*B*B*Y ( 1 )*(F3+ 
2F4)+A*A*Y(3)*SIGEX+B*B*Y(3)*SIGWY 
YDOT ( 4 ) = Y ( 9 ) 

YDOT ( 9 > =-4.*H*H/Z*B**4*< Y ( 4 > -WOT ) +A*A*B*B* ( Y < 3 ) * ( FI 1+F 1 3 ) +8 . *Y I 2 ) * 

1F22 > + 2.*B*C*B*C*Y ( 1 ) *F2+4. *B*B*Y ( 4 ) *S l GWY 
YDOT ( 5 ) =Y ( 10 ) 

YDOT( 1 0 ) = S I GWY/R 

P( 1 )=T 

P< 5)=SIGEX 

RETURN 

END 
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SAMPLE RESULTS (5 MODE MODEL) 
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FOUR MODE MODEL 


// FOR 

* IOCS (1132 PRINTER. CARD. DISK ) 

* LIST SOURCE PROGRAM 
ONE WORD INTEGERS 

FOUR MODE (TWO WAVENUMBER) RUNGE KUTTA SCHfcME 


*********** 


»«***»*•*« 


************* 


A »B 

ITN 

ISHEL 

L 

R 

H 

WTO. WOO. WOT 

Y( 1 ) . . .Y( 5) 

Y (6) . . .Y( 10 ) 
T 

TR 

VO 

DISP 

SIGEX 

PAX 

PCR 

AS.BS.DEL 


AXIAL + CIRCUMFERENTIAL WAVENUMBERS < =M.N ) 

NO. OF WAVENUMBER COMBINATIONS PROVIDED 

SHELL NUMBER 

SHELL LENGTH 

SHELL RADIUS 

SHELL THICKNESS 

IMPERFECTION AMPLITUDES (W20.W11.W02 MODES 
RESPECTIVELY) 

RADIAL DEFLECTIONS 

TIME DERIVATIVES OF Y ( 1 ) ETC. 

TIME PARAMETER ( T I ME /ONE-D I M. WAVESPEED) 
REAL TIME 

AXIAL VELOCITY PARAMETER 
( VELOCITY/ONE-DIM.WAVESPEED) 

AXIAL DISPLACEMENT OF SHELL END 
AXIAL STRESS 
AXIAL LOAD 

CLASSICAL AXIAL BUCKLING LOAD 
STARTING POINT. END POINT AND INCREMENT 
FOR RUNGE-KUTTA CYCLE 


EXTERNAL MYST 
DIMENSION Y ( 10 ) .P ( 5 ) 

RE AD ( 2 • 200 ) ITN 
DO 2 1 = 1. ITN 
READ (2.201) A.B 

200 FORMAT (110) 

201 FORMAT (2F10.5) 

ISHEL = 1 
H=0.0172 

R = 3 . 8 7 A 
WTO=0 • 0 
W00=0. 00001 
WOT = 0.0 
V0=0.0001 
DO 1 J = 1 , 5 
STD = 0 • 0 
P ( 2 ) = VO 
P ( 3 ) = A 
P( A) = 8 

WRITE (3.100) ISHEL . VO. WTO . WOO .WOT . A .B 
100 FORMAT (10H1SHELL NO. .13. /.AH VO.F7.5.&H W20= . F8 . b . 5H Wli=.F8.b. 
15H W02=»F8.5.3H M=,F6.3.3H N= . F6 . 3 . / / . 5X . AHT I ME . 8X . bHWO ( T ) . 6X . 6HW2 
20 (T) .6X.6HW11 (T) .6X.6HW02 (T ) .7X.AHPIT) .5X.8HP (T ) /PCL./ ) 
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Y ( 1 ) =0.0 
Y(2)=WTO 
Y ( 3 ) =WOO 
Y(4)=W0T 
Y ( 5 ) =0.0 

Y ( 6 ) = 0 • 0 

Y ( 7 ) = 0 . 0 

Y ( 8 ) =0.0 

Y ( 9 ) = 0 . 0 
Y ( 10 ) =0.0 
S I GEX = 0 .0 
AS = 0 . 0 

BS = 7 . 0 

DEL = 1 .75 

DO 6 .< = 1.1000 

CALL RKGILI Y . 1 0 . AS . BS . DEL .MYST .P ) 

T = P( 1 ) 

S I GEX = P ( 5 ) 

PAX=2 477000. *R*H*S I GEX 
PCR=1555000.*H*H 
PND=PAX/PCR 
TR=T /70000. 

WRITE(3.101)TR.Y(5) »Y(2) * Y ( 3 ) » Y ( 4 ) .PAX. PND 
101 FORMAT (1X.5IF9.5.3X) .F9.3.3X.F9.5) 

IF (STD-PAX-20. ) 5.1.1 

5 AS = BS 
STD=PAX 

CALL DAT SW (l.KK) 

GO TO (1.6) .« 

6 BS=AS + 7 .0 

1 VO=VO+0. 0001000001 

2 CONTINUE 
CALL EXIT 
FND 

FEATURES SUPPORTED 
ONE WORD INTEGERS 
IOCS 

CORE REQUIREMENTS FOR 

COMMON 0 VARIABLES 82 PROGRAM 456 

END OF COMPILATION 
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n n n n n n 


SUBROUTINE MYSTIM.T.Y.YDOT.P) 

SUBOUTINE MYST FOUR NODE MODEL 

F1.F2... STRESS FUNCTION COEFFICIENTS 

MU POISSON'S RATIO 

SIGEX AXIAL STRESS 

SIGWY CIRCUMFERENTIAL STRESS 

DIMENSION Y< 10 ) i YDOT ( 10) »P ( 5 ) 

VO=P ( 2 ) 

A = P ( 3 ) 

R = P< A ) 

H-0.0172 
R=3.874 
WTO=0 .0 
W00=0, 00001 
WOT = 0 • 0 
NU = 0. AO 

Z = 3.0*( 1.0-NU**2 ) 

L«11.0 

F20=-< 4./R»< Y( 2 )-WTO)+.5*8«B*< Y< 3 )**2-WOO**2 ) ) / ( 16.*A*A> 

Fll=- < A*A/R* ( Y ( 3 ) “WOO ) +2 .*A*A*B*B* i Y ( 2 >*Y ( 3 ) -WTO*WOO+Y ( 3 ) *Y ( 4 l-WOO 
1*W0T ) ) / <A*A+8*B)**2 
F02=-A*A*< Y( 3 )**2-WOO**2 ) / < 32.*8*B) 

F31 = -2.*A*A*B*R*< Y< 2 )*Y< 3 )-WT0*W00) / ( 9 . * A*A+B*B ) **2 
F13=-2 .*A*A*B»B* ( Y ( 3 ) *Y< 4 )-WOO*WOT I / ( A*A+9.*B*B ) **2 
F22=-A*A*B*B* ( Y l 2 ) *Y ( 4 )— WTO*WOT ) / ( A*A+B*B ) **2 
D I SP= VO*T 

SIGEX=NU*SIGWY+DISP/L-A*A*< Y ( 2 ) **2-WTO**2+ . 1 2 5* ( Y < 3 ) **2-WOO**2 ) ) 
SIGWY=NU*SIGEX-B*B*< Y(4)**2-WOT**2+.125*(Y(3)**2-WOO**2) >-Y(5)/R 
YDOT ( 1 1=0.0 
YDOTI 6 )=0.0 
YDOT ( 2 ) = Y ( 7 ) 

YD0T(7)=-4.*H*H/Z*A**4*( Y < 2 ) -WTO > +4 . *A*A/R*F20+A*A*B*B* ( Y(3)*<Fil+ 
1F31 )+8.*Y ( 4 ) *F 2 2 ) +4 . * A*A»Y ( 2 ) *S I GEX 
YDOT ( 3 ) = Y ( 8 ) 

YDOT ( B ) =-H*H/ I 4.*Z ) *( A*A+B*B ) **2* ( Y t 3 ) -WOO ) +A* A / R*F 1 1 +2 • *A*A*6*B* ( 
1Y(2>*<F11+F31>+Y<3)*(F20+F02H-Y(4)*<FI1+F13) ) +A*A*Y t 3 > *S IGEX+B*B*Y 
2(3) *S I GWY 
YDOT ( 4 ) = Y ( 9 ) 

YD0T(9)=-4.*H«H/Z*B**4*< Y ( 4 > -WOT ) +A*A*B*B* ( Y ( 3 ) * ( F 1 1+F 1 3 ) +B . *Y ( 2 ) * 
1F22)+4.*R*B*Y(4)*SI GWY 
YDOT ( 5 >=Y( 10) 

YDOT ( 1 0 ) =S I GW Y/R 

P( 1)=T 

P ( 5 ) =S IGEX 

RETURN 

END 
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SAMPLE RESULTS (4 MODE MODEL) 


SHELL NO . 1 


V 0 = 0. 00010 

W 20 = 0.00000 

Wll = 0.00001 

W 02 = 0.00000 M > 3.410 

N * 3.410 


TIME 

worn 

W 20 IT ) 

W 1 1 < T ) 

W 02 IT ) 

PIT ) 

P ( TI/PCL 

0.00010 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

10.303 

0.02283 

0.00020 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

21.006 

0.04566 

0.00030 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

31.509 

0.06849 

0.00040 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

42.012 

0.09132 

0.00050 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

52.515 

0.11415 

0.00060 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

63.018 

0.13698 

0.00070 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

73.521 

0.15981 

0.00080 

- 0.00000 

- 0,00000 

0.00001 

- 0.00000 

84,023 

0.18265 

0.00090 

- 0,00000 

- 0.00000 

0.00001 

- 0.00000 

94.528 

0.20548 

0.00100 

- 0,00000 

- 0.00000 

0,00001 

- 0,00000 

105.031 

0.22831 

0.00110 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

115.534 

0.25114 

0.00120 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

126.037 

0.27397 

0.00130 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

136.540 

0.29680 

0.00140 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

147.043 

0.31963 

0.00150 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

157.547 

0.34247 

0.00160 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

168.050 

0.36530 

0.00170 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

178.553 

0.38813 

0.00180 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

189.056 

0.41096 

0.00190 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

199.559 

0.43379 

0.00200 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

210.062 

0.45662 

0.00210 

- 0.00000 

- 0.00000 

0.00002 

- 0.00000 

220.565 

0.47945 

0.00220 

- 0.00000 

- 0.00000 

0.00002 

- 0.00000 

231.069 

0.50229 

0.00230 

- 0.00000 

- 0.00000 

0.00002 

- 0.00000 

241.572 

0.52512 

0.00240 

- 0.00000 

- 0.00000 

0.00002 

- 0.00000 

252.075 

0.54795 

0.00250 

- 0.00000 

- 0.00000 

0.00002 

- 0.00000 

262.578 

0.57078 

0.00260 

- 0.00000 

- 0.00000 

0.00002 

- 0.00000 

273.081 

0.59361 

0.00270 

- 0.00000 

- 0.00000 

0.00002 

- 0.00000 

283.584 

0.61644 

0.00280 

- 0.00000 

- 0.00000 

0.00003 

- 0.00000 

294.087 

0.63927 

0.00290 

- 0.00000 

- 0.00000 

0.00003 

- 0.00000 

304.590 

0.66210 

0.00300 

- 0.00000 

- 0.00000 

0.00003 

- 0.00000 

315.093 

0.68494 

0.00310 

- 0.00000 

- 0.00000 

0.00004 

- 0.00000 

325.596 

0.70777 

0.00320 

- 0.00000 

- 0.00000 

0.00004 

- 0.00000 

336.099 

0.73060 

0.00330 

- 0.00000 

- 0.00000 

0.00005 

- 0.00000 

346.603 

0.75343 

0.00340 

- 0,00000 

- 0.00000 

0.00005 

- 0.00000 

357.106 

0.77626 

0.00350 

- 0.00000 

- 0.00000 

0.00006 

- 0.00000 

367.609 

0.79909 

0.00360 

- 0.00000 

- 0.00000 

0.00007 

- 0.00000 

378,111 

0.82192 

0.00370 

- 0.00000 

- 0.00000 

0.00008 

- 0.00000 

388.614 

0.84475 

0.00380 

- 0.00000 

- 0.00000 

0.00010 

- 0.00000 

399.116 

0.86758 

0.00390 

- 0.00000 

- 0.00000 

0.00012 

- 0.00000 

409.618 

0.89041 

0.00400 

- 0.00000 

- 0.00000 

0.00015 

- 0.00000 

420.120 

0.91324 

0.00410 

- 0.00000 

- 0.00000 

0.00019 

- 0.00000 

430.619 

0.93606 

0.00420 

- 0.00000 

- 0.00000 

0.00025 

- 0.00000 

441.116 

0.95888 

0.00430 

- 0.00000 

- 0.00001 

0.00034 

- 0.00000 

451.606 

0.98168 

0.00440 

- 0.00000 

- 0.00004 

0.00048 

- 0.00000 

462.078 

1.00445 

0.00450 

- 0.00000 

- 0.00010 

0.00070 

- 0.00000 

472.501 

1.02710 

0.00460 

- 0.00000 

- 0.00027 

0.00107 

- 0.00001 

482.723 

1.04932 

0.00470 

- 0.00001 

- 0.00075 

0.00175 

- 0.00002 

491.821 

1.06910 

0.00480 

- 0.00003 

- 0.00216 

0.00321 

- 0.00007 

492.640 

1.07088 

0.00490 

- 0.00015 

- 0.00606 

0.00718 

- 0.00022 

431.792 

0.93861 
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shell no . i 

V 0 * 0. 00020 W 20 = 0.00000 Wll = 0.00001 W 02 = 0.00000 M = 3.410 N = 3.410 


TIME 

WO ( T ) 

W 20 ( T ) 

Wll (T ) 

W 02 ( T ) 

PIT ) 

PIT) /PCL 

0.00010 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

21.006 

0.04566 

0.00020 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

42.012 

0.09132 

0.00030 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

63.018 

0.13698 

0.00040 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

84.025 

0.18265 

0.00050 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

105.031 

0.22831 

0.00060 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

126.037 

0.27397 

0.00070 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

147,044 

0.31963 

0.00080 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

168.050 

0.36530 

0.00090 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

189.056 

0.41096 

0.00100 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

210.062 

0.45662 

0.00110 

- 0.00000 

- 0.00000 

0.00002 

- 0.00000 

231.069 

0.50229 

0.00120 

- 0.00000 

- 0.00000 

0.00002 

- 0.00000 

252.075 

0.54795 

0.00130 

- 0.00000 

- 0.00000 

0.00002 

- 0.00000 

273.081 

0.59361 

0.00140 

- 0.00000 

- 0.00000 

0.00002 

- 0.00000 

294.087 

0.63927 

0.00150 

- 0.00000 

- 0.00000 

0.00003 

- 0.00000 

315.094 

0.68494 

0.00160 

- 0.00000 

- 0.00000 

0.00004 

- 0.00000 

336.100 

0.73060 

0.00170 

- 0.00000 

- 0.00000 

0.00004 

- 0.00000 

357.106 

0.77626 

0.00180 

- 0.00000 

- 0.00000 

0.00006 

- 0.00000 

378.112 

0.82192 

0.00190 

- 0.00000 

- 0.00000 

0.00007 

- 0.00000 

399.118 

0.86758 

0.00200 

- 0.00000 

- 0.00000 

0.00010 

- 0.00000 

420.123 

0.91325 

0.00210 

- 0.00000 

- 0.00000 

0.00013 

- 0.00000 

441.127 

0.95890 

0.00220 

- 0.00000 

- 0.00000 

0.00019 

- 0.00000 

462.128 

1.00456 

0.00230 

- 0.00000 

- 0.00001 

0.00030 

- 0.00000 

483.122 

1.05019 

0.00240 

- 0.00000 

- 0.00003 

0.00047 

- 0.00000 

504.093 

1.09578 

0.00250 

- 0.00000 

- 0.00011 

0.00081 

- 0.00000 

524,974 

1.14117 

0.00260 

- 0.00000 

- 0.00041 

0.00148 

- 0.00001 

545.304 

1.18536 

0.00270 

- 0.00003 

- 0.00167 

0.00300 

- 0.00006 

559.575 

1.21638 

0.00280 

- 0.00014 

- 0.00688 

0.00746 

- 0.00021 

483.799 

1.05166 


SHELL NO . 
V 0 = 0. 00030 

1 

W 20 = 0.00000 

Wll = 0.00001 

W 02 « 0.00000 M » 3.410 

N “ 3.410 


TIME 

WO I T ) 

W 20 I T ) 

Wll ( T ) 

W 02 I T ) 

PIT ) 

PIT) /PCL 

0.00010 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

31.509 

0.06849 

0.00020 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

63.018 

0.13698 

0.00030 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

94.528 

0.20548 

0.00040 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

126.037 

0.27397 

0.00050 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

157.547 

0.34247 

0.00060 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

189.056 

0.41096 

0.00070 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

220.566 

0.47945 

0.00080 

- 0.00000 

- 0.00000 

0.00002 

- 0.00000 

252.075 

0.54795 

0.00090 

- 0.00000 

- 0.00000 

0.00002 

- 0.00000 

283.584 

0.61644 

0.00100 

- 0.00000 

- 0.00000 

0.00003 

- 0.00000 

315.094 

0.68494 

0.00110 

- 0.00000 

- 0.00000 

0.00004 

- 0.00000 

346.603 

0.75343 

0.00120 

- 0.00000 

- 0.00000 

0.00005 

- 0.00000 

378.112 

0.82192 

0.00130 

- 0.00000 

- 0.00000 

0.00007 

- 0.00000 

409.621 

0.89042 

0.00140 

- 0.00000 

- 0.00000 

0.00010 

- 0.00000 

441.129 

0.95891 

0.00150 

- 0.00000 

- 0.00000 

0.00015 

- 0.00000 

472.636 

1.02740 

0.00160 

- 0.00000 

- 0.00000 

0.00024 

- 0.00000 

504.136 

1.09587 

0.00170 

- 0.00000 

- 0.00001 

0.00041 

- 0.00000 

535.617 

1.16430 

0.00180 

- 0.00000 

- 0.00008 

0.00077 

- 0.00000 

567.010 

1.23254 

0.00190 

- 0.00000 

- 0.00039 

0.00158 

- 0.00001 

597.780 

1.29943 

0.00200 

- 0.00003 

- 0.00207 

0.00363 

- 0.00007 

618.704 

1.34491 

0.00210 

- 0.00024 

- 0.01052 

0.01051 

- 0.00028 

423.683 

0.92098 
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SHELL NO . 
V0 = 0. 00050 

1 

W20 = 0.00000 

Wll = 0.00001 

W02 « 0.00000 

M = 3.410 IM = 

3.410 

TIME 

worn 

W20 ( T ) 

W 1 1 ( T ) 

W02 ( T ) 

PIT ) 

0 .00010 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

52.513 

0.00020 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

105.031 

0.00030 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

137.347 

0,00040 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

210.062 

0.00050 

- 0.00000 

- 0.00000 

0.00002 

- 0.00000 

262 ,578 

0.00060 

- 0.00000 

- 0.00000 

0. 00002 

- 0.00000 

315.094 

0.00070 

- 0.00000 

- 0.00000 

0.00003 

- 0.00000 

367.609 

0.00080 

- 0.00000 

- 0.00000 

0.00003 

- 0.00000 

420.123 

0.00090 

- 0.00000 

- 0.00000 

0.00008 

- 0.00000 

472.639 

0.00100 

- 0.00000 

- 0.00000 

0.00014 

- 0.00000 

525.152 

0.00110 

- 0.00000 

- 0.00000 

0.00028 

- 0.00000 

577.634 

0.00120 

- 0.00000 

- 0.00003 

0.00058 

- 0.00000 

630 .103 

0.00130 

- 0.00000 

- 0.00021 

0.00137 

- 0.00001 

682 .137 

0.00140 

- 0.00003 

- 0.00167 

0.00365 

- 0.00006 

726.332 

0 . 0-0150 

- 0.00028 

- 0.01277 

0.01232 

- 0.00030 

439.403 


SHELL NO . 
VOsO . 00040 

1 

W 20 = 0.00000 

Wll = 0.00001 

W 02 = 0.00000 

M » 3.410 N = 

3.410 

TIME 

WO ( T ) 

W 201 T ) 

WU IT ) 

W 02 ( T ) 

PIT ) 

0.00010 

- 0.00000 

- 0.00000 

0.00001 

- 0.00000 

42.012 
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- 0.00000 
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0.00001 
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602.650 


PITl/PCL 

0.11415 

0.22831 

0.34247 

0.43862 

0.8707b 

0.66494 

0.79909 

0.91328 

1102740 

1.14133 

1.2336b 

1.36969 

1.46263 

1.37931 

0.95316 


PITl/PCL 

0.09132 

0.18265 

0.27397 

0.36530 

0.43662 

0.54795 

0.63927 

0.73060 

0.62192 

0.91323 

1.00457 

1.09589 

1.18719 

1.27640 

1.36906 

1.45271 

1.31002 
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APPENDIX C: DYNAMIC PROPERTIES OF EPOXY PLASTIC 

In Chapter 3 of the main text, it was explained that because of the 
particular requirements of the spin-casting technique, the shell specimens were 
manufactured from an epoxy plastic material. However, a well recognized feature 
of polymer solids such as epoxies is that their stress-strain characteristics 
are normally rate dependent, and that in extreme cases, their response to dynami- 
cally applied loads can differ very signficantly from their static behaviour. 

The visco-elastic properties of the epoxy mixture used in the manufacture 
of the shell specimens were determined in a series of tests which employed three 
experimental approaches. The first and simplest technique involved the use of a 
small pneumatic testing rig in which several small samples of the material were 
tested in tension and compression at strain rates varying from the quasi-static 
up to 1,5 x 10"^ sec -1 . In the second series of experiments, considerably higher 
strain rates were achieved with a split Hopkinson bar apparatus, while by the 
third method, material properties were determined from the transmission charac- 
teristics of high-frequency stress waves initiated by an explosion. The first 
two sets of experiments were performed by the author, while the third was con- 
ducted by Zimeik^ 2 , 

The techniques and principles employed in the first series of tests 
were quite conventional and straightforward. The loading elements of the. test- 
ing machine consisted of a pair of pneumatic cylinders so that the loading rate 
of the machine corresponded to the rate at which the cylinders were pressurized. 
Strain was measured by small foil resistance strain gauges bonded directly to 
the sample while the load was measured by a simple strain gauge load cell 
mounted in series with the sample piece between the platens of the testing 
machine. The strain gauges on the sample and the load cell were connected to 
bridge circuits, the output signals from which were displayed on a dual beam 
oscilloscope. The results obtained with this arrangement are shown in Figure 31 • 

The split Hopkinson bar technique used to measure material properties 
at medium and high strain ra.tes was somewhat more complicated in concept than the 
tests described above. The apparatus consisted essentially of two slender bars 
(15 inches long, ,375 inches in diameter) mounted in tandem on a common axis 
with a small (,5 inches long by ,375 inches in diameter) sample of the material 
sandwiched between them (see Figure 30-) , 

In operation, the free end of one of the bars (termed the weigh bar) 
was struck by a projectile fired from a large air gun so that a compressive one- 
dimensional stress wave propagated down the bar and through the sample into the 
downstream or anvil bar. Foil resistance-type strain gauges were mounted on 
each bar close to the sample interfaces to measure the incoming, reflected and 
transmitted stress waves. The strain signals were displayed on an oscilloscope 
and recorded photographically. By assuming that the stress waves were essentially 
one-dimens icnal (due to the geometry of the bars) it was possible to calculate . 
the particle velocities in the two bars from the waveforms and hence to derive 
the strain rate and strain from the relative velocities and displacements of 
the two bar-sample interfaces . The sample stress was taken as being equal to 
the average of the stresses in the two bars. 

The mathematical model used to calculate the strains was developed 
as follows. An element of bar of length Ax is shown below. The axial stresses 
and displacements are represented by cr(x,t) and u(x,t) respectively, while the 
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mass density of the bar material's ^epres^nt^ ^ p. 

r* r* 



The unbalanced force on the element will be equal to the mass times 
acceleration of the element, so that in the limit, as Ax approaches zero we 
may write 


dcr(x,t) _ pd u(x,t) 

" ' St 2 

If it is assumed that the material is linearly elastic, that is 
cr(x,t) = Ee(x,t) = 


C.l 


C .2 


then C.l becomes 


„ ^2 -\2 

E o u o u _ n 

p a x 2 at 2 


c.3 


This is a wave equation with the solution 


where c 



u(x,t) = f(x + ct) + g(x-ct) 


C.k 


If we consider the case of a stress wave of magnitude 2 moving to the right along 
a previously undisturbed bar, we have 

cr(x,t) = E g' (x-ct) = 0 C.5 

in the undisturbed region and 

cr(x,t) = £(x,t) = E g’(x-ct) C.6 

behind the wave front. The particle velocity behind the wave front is 

v( X ,t) = ^ (x,t) = - c g'(x-ct) 

so that 

2(x,t) = pcv(x,t) C.7 
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When the stress wave encounters an interface at X between two different 
materials, a more complicated situation arises with displacements due to the 
incident wave g (x-ct), a reflected wave f (x + ct) and a wave transmitted into 

xA. 

the new material Ug(x,t). At the interfaces, continuity demands that the dis- 
placements and velocities in two materials be equal, that is, 


g A (X-ct) + f A (X + ct) = Ug(X,t) 


-cg’(X-ct) + cf* (X + ct) = 


dU B 

'ST 


(X,t) 


c.8 

c.9 


Also, the equilibrium condition specifies that the stresses in both materials 
be equal: 


E (g^ (X-ct) + f' A (X-ct)) = cr B (x,t) 


C.10 


An important design feature of the split Hopkins on bar apparatus is 
that the incident stress wave is of essentially constant magnitude for the 
duration of the pulse. If this stress level is written £ , then 

Z 

g^ (X-ct) = ^ C.ll 


If this stress level is known and' the total stress level at the interface is 
determined from strain gauge measurements, we can write 

a (X,t) 

q(x-ct)=^ £ a €.12 

The stress at the interface cr is, of course, equal to cr 

1 £> 

Thus, the velocity in material B adjacent to the interface is 


du (X,t) 

W ” E (X,t) ' 2Z A ) 


(ct (X,t) - 2£ a ) 
pc 


c.13 


At the sample-anvil bar interface, waves transmitted from the sample 
into the bar are entering an undisturbed medium so that it is possible to use 
Equation C.7 5 written here as 


dt 


(x + £,t) 


a (x + i,t) 

pc 


c.i4 


The sample strain rate is then 
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(X + i,t) 


C . 15 


£(t) 


1 


l 


r au B 

5t~ 


(x 3 t) 


5u b 


and the sample strain at time T is 


£(T) = 


T 


£ ( t ) dt 


C.l6 


It should be noted that these equations make no assumptions as to the 
material properties of the sample itself and are applicable regardless of any 
energy dissipation or internal reflections within the sample piece . For prac- 
tical reasons, it was impossible to measure the bar strains' precisely at the 
interfaces as was implied in the above derivation. If, however, it is assumed 
that there is negligible dispersion or attenuation in the steel bars between 
the ends and the strain gauge stations (a distance of approximately 3/8 of an 
inch) , the equations remain completely applicable,, 

An attractive feature of the split Hopkinson bar technique is that 
the instrumentation is carried entirely on the apparatus, so that the test 
specimens are very simple and economical to produce „ Nevertheless, several 
sample pieces were specially cast with small strain gauges embedded directly 
at their centre. Comparison between the internal strains measured directly 
by this technique and the strains calculated from the wave data indicated that 
the one-dimensional wave analysis outlined above was capable of producing quite 
accurate results. 

A full description of the third experimental method used in determin- 
ing dynamic stress-strain properties of the epoxy material is contained in 
Reference 32. Briefly, this method was based on an analysis of the wave 
transmission properties of the material. A long slender sample was mounted 
in linear bearings, which permitted axial motion, and loaded impulsively at 
one end by a small explosive charge. Strain gauges mounted at the mid-point 
of the bar were used to measure the strain pulse as it travelled back and 
forth along the bar and was reflected at the free ends. Each of these strain 
pulses was represented mathematically as a Fourier series so that by comparing 
the magnitude of the Fourier components of each pulse, it was possible to de- 
termine the attenuation coefficient and phase velocity of a strain wave as a 
function of frequency. This data was used to determine the constitutive prop- 
erties of the material in terms of a complex Young’s Modulus (E^+ iE^) where 
E^ is the elastic stiffness and E^ is representative of the energy dissipation 

in. the material. 

These three experiments were complemented by earlier work done by 
Tennyson 11 , in which measurements of the static modulus and Poisson's ratio 
were made. 


Sample results from the three dynamic experiments are plotted to- 
gether in Figure 32, while the results from the first low strain rate series 
are shown in Figure 31° The results from the wave transmission experiment 
are plotted as a straight line. This must not be taken to imply that this 
experiment predicts completely linear behaviour at all stress levels, but 
rather it simply reflects the fact that at the strain levels encountered here 
and because of the nature of the analysis procedure, no information on yield 


C4 



points or non-linear behaviour was obtained. 

During the course of the dynamic shell buckling tests, the maximum 
pre-buckling strains were of the order of 4 x 10 “ 3 , while the maximum pre- 
buckling strain rates were approximately 2.5 sec - '*- and the maximum buckling 
strain rates equal to about 4 sec - *-. This range of strains corresponds to a 
region of the stress- strain diagrams where the results obtained with all three 
experimental methods are in good agreement. It was decided, therefore, that 
the material could be considered as being linearly elastic with a Young's modulus 
of 4.05 x 10 ^ psi over the entire range of strain rates encountered in the shell 
buckling tests. It was also assumed that the statically determined Poisson's 
ratio of 0.4 would remain constant under these conditions. 


C5 



SERIES I: NEAR-PERFECT SHELLS 
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TABLE 2 


COMPARISON OF THEORETICAL AND EXPERIMENTAL BUCKLING WAVELENGTHS 


SHELL 3 (Figure 28) 

AXISYMMETRIC MODE (W(t) cos2Mx) 

Theoretical Wavelength 
(inches) 

-924 

ASYMMETRIC MODE (W(t) cosMxcosNy) 

Theoretical Wavelengths 
(inches) 


Experimental Wavelength 
(inches) 

.93 


Experimental Wavelength 
(inches) 


Axial Circumferential 


Axial Circumferential 


1.85 1.62 


2.0 2.1 


SHELL 9 (Figure 29) 

AXISYMMETRIC MODEL (W(t) cosKx) 

Theoretical Wavelength 
(inches) 

2.756 

ASYMMETRIC MODE (W(t)cosMxcosNy) 

Theoretical Wavelengths 
(inches) 

Axial Circumferential 

1.96 1.96 


Experimental Wavelength 
(inches) 

2.3 


Experimental Wavelengths 
(inches) 

Axial Circumferential 

1.7 1.3 





Quasi -Random Imperfection 


FIG. 1(b) IMPERFECTION GEOMETRIES 


FIG. 2 COMPUTER GENERATED RADIAL DISPLACEMENTS 8 AXIAL LOADS 
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FIG. 3 
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FIG. 5 PHOTOGRAPH OF PROFILE MEASURING APPARATUS 





FIG. 6 TESTING MACHINE AND INSTRUMENTATION 
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G. 7 SCHEMATIC VIEW OF TESTING MACHINE 
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FIG. 9 OSCILLOSCOPE PHOTOGRAPHS FROM A TYPICAL 
TEST SERIES (SHELL 8) 
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FIG. 10(a) THEORETICAL AXIAL LOADS VS TIME 
PERFECT SHELL MODEL 
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FIG. 10(b) THEORETICAL BUCKLING LOADS VS 
DEFLECTION RATE 



VARIATION IN DYNAMIC BUCKLING LOADS WITH IMPERFECTION AMPLITUDE 



FIG. 12 VARIATION IN DYNAMIC BUCKLING LOADS WITH IMPERFECTION 
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FIG. 14 


SHELL MEDIAN SURFACE PROFILES 



FIG. 15 SHELL MEDIAN SURFACE PROFILES 
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Shell 6 



FIG. 16 SHELL WALL a THICKNESS PROFILES 




G. 17 SHELL WALL & THICKNESS PROFILES 




DYNAMIC BUCKUNG LOADS VS AXIAL DEFLECTION RATE 
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FIG. 20 DYNAMIC BUCKLING LOADS VS AXIAL DEFLECTION RATE 
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FIG. 21 DYNAMIC BUCKLING LOADS VS AXIAL DEFLECTION RATE 
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FIG. 22 DYNAMIC BUCKLING LOADS VS AXIAL DEFLECTION RATE 
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FIG. 23 DYNAMIC BUCKLING LOADS VS AXIAL DEFLECTION RATE 
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FIG. 24 DYNAMIC BUCKLING LOADS VS AXIAL DEFLECTION RATE 
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FIG. 25 DYNAMIC BUCKLING LOADS VS AXIAL DEFLECTION RATE 



FIG. 26 DYNAMIC BUCKLING LOADS VS AXIAL DEFLECTION RATE 



Fourier Component Approximation 
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SHELL 10, QUASI - RANDOM IMPERFECTION 
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FIG. 29 HIGH-SPEED PHOTOGRAPHS OF DYNAMIC BUCKLING 
PROCESS 

Shell 9 Vo = 2 x I0" 4 
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FIG. 33 STRESS-STRAIN CURVES FOR EPOXY 


PLASTIC AT HIGH STRAIN RATES 



STRAIN 


